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Summary 


The theory of the propagation of slowly varying electromagnetic 
disturbances through partially ionized gas is developed and applied to 
the Earth’s atmosphere and interplanetary space. The medium must 
be regarded as two separate, co-existing gases, an electron-ion plasma 
and neutral atoms which move to some extent independently. Quantita- 
tive results are given for a model atmosphere out to several Earth radii: 

(1) Up toa few hundred kilometres the medium behaves, for waves 
of all periods between 1s and 10's, as a rigid conductor and as a disper- 
sive medium. 

(2) Above about 10km the disturbances travel as hydromagnetic 
waves in the ion plasma alone. Losses are small in this region but 
transmission is likely to be complicated by anisotropic transmission of 
the O wave and by refraction and partial reflection. 

(3) The currents responsible for all observed geomagnetic dis- 
turbances must flow at levels below about 1 000 km. 

(4) The problem of the penetration of solar gas into the Earth’s 
field is discussed. 

(5) Some earlier theories of the main phase of an SC storm are 
discussed and an alternative suggested, whereby some lines of force 
of the Earth’s field are carried away along the Sun-Earth line. 

(6) Some properties of micropulsations are explained in terms of 
the transmission of hydromagnetic waves. 

(7) The “effective” conductivity (o3) of the atmosphere out to 
several Earth radii has its maximum about 100 km, above which it rapidly 
falls to about 10~!%e.m.u. or less and does not rise again. 


1. Introduction 


Geomagnetic disturbances observed near the Earth’s surface have fluctuation 
periods ranging from a fraction of a second in the case of some micropulsations 
up to some days for the slowest storm variations. 
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It is generally accepted that some of these disturbances originate outside the 
Earth’s atmosphere. They must comprise, then, electromagnetic waves which 
propagate through the Earth’s atmosphere and perhaps through interplanetary 
space. In so propagating they will suffer phase delay, refraction, absorption and 
partial reflection and the disturbance finally observed at the Earth’s surface 
will depend on all these factors as well as on the nature of the original 
disturbance. 

Propagation through the ionosphere has been considered by Ashour & Price 
(1948) and others in a theory in which the medium does not move under the 
influence of the electromagnetic field. A number of workers have considered the 
passage of micropulsations as downward travelling hydromagnetic waves. As 
will be seen below, these theories do not take full account of the existence of two 
gas components, an ion-electron plasma (here called the plasma) and neutral 
atoms. These gases together constitute the whole, partially ionized, gas of the 
atmosphere and interplanetary space. Under the influence of an electromagnetic 
field they move to some extent, independently. 

The transmission of geomagnetic disturbances is a problem in hydromagnetics 
in a gas mixture which may not be treated as a single fluid. It is found that the 
concept of electrical conductivity of the whole gas is not valid and that the medium 
is dispersive for hydromagnetic waves. 

The theory developed here must be regarded as a preliminary step. It con- 
cerns weak waves in a uniform medium and so excludes possibly important 
refraction and reflection effects. Nevertheless, it may reveal some significant 
transmission effects and throw some light on the difficult theory of the origins of 
the geomagnetic disturbances. 


2. The transmission equations 


The transmission of hydromagnetic waves through a partially ionized gas may 
be studied by using Maxwell’s field equations, together with the separate equations 
of motion of the three gas components present (for application to micropulsations 
and references to earlier work, see Akasofu 1956). This approach is of course 
valid, but adds greatly and unnecessarily to the complexity of the problem. The 
reason is that Maxwell’s displacement currents may be neglected in these problems 
and so, to the same order of accuracy, we may write divj = 0, where j is the 
electric current density.* This means that the electron gas and the heavy ion gas 
may be considered a single dynamic entity and have a single force-momentum 
equation (for a previous discussion see Piddington 1954). 

The situation is much simplified by using the equations of motion of two gases 
only: the plasma and the neutral atom gas. It is further simplified by neglecting 
ordinary viscosity due to thermal motions of the particles; this is usually unim- 
portant compared with ‘“‘neutral atom damping” (Piddington 1957). Finally, the 
effects of gas pressure gradients caused by the wave are neglected. These are not 
present when the waves are propagated along the magne -c field; for transverse 
propagation compression effects are small provided the velocity of sound is small 
compared with the hydromagnetic velocity V (Piddington 1955), which is the case, 
except for waves of long period in ionospheric E region. 


* This does not mean that electric charge may not accumulate on boundaries, but that the 
currents involved in causing such accumulations have negligible dynamic effects. 
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The hydromagnetic equations may now be listed (Piddington 1956): 
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Here 6}, 62, 63 are the well-known conductivity components (e.m.u.) of the plasma 
alone and H the magnetic field, with a steady component Ho directed along O.z 
The velocity vectors of the plasma and neutral atom gas are v and v’ and 7 = pp/p 
is the ratio of the mass densities of neutral atoms and ions. Finally 7 is the period 
taken to accelerate the neutral atoms to the same velocity as the plasma, assuming 
the latter remains constant during this period (see Piddington 1957, Appendix I). 

Although we are dealing with a partially ionized gas, equation (1) refers to a 
fully ionized gas because the effects of neutral atoms are introduced separately by 
equations (2) and (3). This method is essential when dealing with waves of fre- 
quency greater than the neutral atom collision frequency, when the concept of 
conductivity in the partially ionized gas is not valid. This applies for most geo- 
magnetic disturbances. 

Equation (1) was derived with the simplifying assumption that the collision 
rate of electrons with heavy ions was substantially greater than with neutral atoms, 
an assumption which is certainly not valid in the lower ionosphere. However, the 
effects of collisions with neutral atoms remain unimportant. The losses are 
increased by such collisions (by a factor of ~ 10% in the E region) but the corres- 
ponding losses due to collisions of heavy ions and neutral atoms are much greater, 
the factor being ~ 10° (compare oo and o3 in Table 1 below). 

Equations (2) and (3) are linearized, being valid only when the perturbations 
v, v and Ho-H are small compared with V and Hp. This condition is satisfied 
generally, except at several Earth radii where the field may be violently disturbed. 


3. The transmission parameters 
Solutions of the above equations are sought, of the form exp i(wt—kx). The 
dispersion equation is (Piddington 1956) 
w2 i 


— = OV2S?+-, 
un? : (5) 


where O = {1+7(1+iw7)}" (6) 


and S has very nearly the values cos @ and unity for the ordinary (O) and extra- 
ordinary (E) waves respectively (Piddington 1955), @ being the angle between the 
wave normal and Hy. The term i/e corresponds to absorption due to the finite 
value of oo (that is, due to electron-ion collisions) and as seen in the Appendix, 
may generally be neglected. When this is done the dispersion equation becomes 


k= 1+7(1+iwr)}}, (7) 
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There are three limiting cases, corresponding to different physical situations: 
(a) Waves of very low frequency, wr > 1. The dispersion equation reduces to 
1+7)* i 
a w (1 | (8) 
VS 2(1+7) 


which describe hydromagnetic waves travelling with phase (and group) velocity 
VoS, where 





Vo = V(1+n)7. (9 
The absorption coefficient « (the wave decays as exp(—«x)) is given by 


9 


k= SE at alle (10) 
2VoS(1 +7) 


The absorption coefficient of a hydromagnetic wave in a single gas of conduc- 
tivity o is (Piddington 1955) 


Kk = w*(8r0S3Vo3)-1. (11) 
Putting S = 1 and using equation (g), the two expressions for « are identical, 
provided 
(p+ pn)? 
c= en, (12) 
pntH2 


which is the value of o3 for the whole gas (see Appendix, the o9~! term being 
neglected). Hence the waves described by equations (9) and (10) travel in the 
composite gas, which may be considered as a fluid of density (p+ pn) and effective 
conductivity o3 as found in the Appendix. 


(b) A second limiting case occurs when 7 > wr > 1. The dispersion equation 
then reduces to 


k ~ a3 (2) ¢-a, (13) 


which describes waves travelling with velocity VoS(2wr) and suffering absorption 


I =) 
k= 73(~ . (14) 


These waves travel neither in the medium as a whole (their velocity would be 
VoS) nor in the plasma alone. Their velocity corresponds to that in a single 
medium of density py,/2w7r, which lies between py and p and varies with the wave 
frequency so that the medium is dispersive. 

Another way of looking at this case is to combine equations (9), (12) and (14) 
from which, remembering that p, > p, we have 


K = (2703w)*S-. (15) 
Apart from the factor S-, this is identical with the solution of the equation 
oH 
V*H = 4703—, (16) 


ct 


for waves of the form required here. This field equation, recognizable from 
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“skin effect” theory, applies when the conductor does not move. Thus equation 


(15) corresponds, apart from geometrical complexities, to the situation assumed 
by Ashour & Price (1948). 


(c) A third case of interest is that for which wr > 1+, when 


w in 
ean F me ) 
VS (: aa (17) 


These waves have velocity VS so that they are travelling in the plasma alone. 
The absorption coefficient is 


i] 
in i 18 
<" 2VSr (18) 


so that the neutral atom gas is important in determining the absorption rate. 


4- A model of the atmosphere and interplanetary medium 


In order to permit a quantitative discussion of the results of the previous 
section a (very incomplete) model of the Earth’s atmosphere* is given in Table 1. 
Included are the number densities N, Ny of the heavy ions and neutral atoms, their 
average atomic weights A (assumed equal) and the gas temperature 7. 

At the higher levels 4,, h,, h,, whose actual heights are unspecified but may be 
in the range 3-8 Earth radii, the true conditions are virtually unknown. Level A, 


may be that through which “whistling atmospherics” pass and so have an electron 
(and ion) density ~ 108cm-%. The temperature is likely to be controlled largely 
by conductivity downwards and so should be between ~ 10°°K in the ionosphere 
and perhaps ~ 105°K in interplanetary space. The corresponding proportion of 
neutral atoms should be ~ 1 in 200 (Nikolsky 1958). 

The two columns headed h,, h, represent alternative possibilities of conditions 
at the border between the Earth’s atmosphere and true interplanetary space. The 
corresponding magnetic field strength listed is that at ~ 7 Earth radii under static 
conditions. The values of N = 10%cm-% is probably an upper limit and there is 
evidence that it may be too large (Blackwell 1957). Consequently a second alterna- 
tive N = 10%cm~% is also listed. The two corresponding values of N, are chosen 
to provide probable upper and lower limits of o3. With the maximum value of N 
is associated the minimum possible value of N, corresponding to a high tempera- 
ture in a pure hydrogen atmosphere. The temperature may be as high as ~ 10°°K 
(Chapman 1957) in which case a state of equilibrium between ionization by 
collision and recombination gives N,/N ~ 10-6 (Shklovsky 1951). The ratio is 
independent of the gas density. The maximum likely value of N, is associated with 
the lower value of N (at h,). The gas temperature may be lower and dust may 
also play a part in increasing N,/N to about 1/500 (Shklovsky 1958). 

Also listed in Table 1 are the values of magnetic pressure, Ho?/87, and gas 
pressure at the three upper levels. In each case the magnetic pressure is sub- 
stantially larger, so that the approximation of neglecting elastic forces in the gas 


compared with those in the field (neglecting the velocity of sound compared 
with V) is a reasonable one. 


* I wish to thank Dr D. F. Martyn and Dr G. R. Ellis for access to data averaged from 
available experimental and theoretical sources. 
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Table 1 


A model atmosphere 


(km) 100 ha 


(cm~%) 10° 
(cm~) 3 X 1012 
30 

(°K) 250 

(s) 2-0 X 104 1-3 X10? 

(gauss) 0°33 . 0-012 

(e.m.u.) — Io X 1078 

(e.m.u.) 6-7 x 10714 1-7 X 10722 

(cm s~) 4°2 X10? 8-3 x 107 6-9 x 108 2°2 X107 
Vo (cm s~!) 7°6 x 108 8-3 x 107 6-9 x 108 2°2 X 107 
H,?2/8m (dyn cm-*) — 5°7 X 1078 4°0 X 1078 4°0 x 1078 
NkT (dyn cm~*) _— 1-4 X 107° 1-4 X1078 >1+4 X 1078 


o~I°3 x 108 
1073 


The table also gives oo, the conductivity along the magnetic lines of force; this 


is the effective transverse conductivity in the absence of all neutral atoms (when 
69 = 53, see Appendix). 


5. The transmission of geomagnetic disturbances 


When the data of Table 1 are substituted in the equations of the previous 
section the propagation constants of electromagnetic waves through the atmo- 
sphere may be determined. This has been done for three wave periods, 1s, 


100s and 100008, which cover most of the range of observed fluctuations. For 
each period and for each level of the model atmosphere the propagation velocity 
(V(E), etc.) is determined and the distance (8(Z) etc.) in which the wave is 


attenuated by a factor e~!. The latter variable is the inverse of «. The results are 
given in Table 2. 


Table 2 


Wave velocities and rates of absorption 


wave period (s) 100 10000 


w 6°3 X 1072 6-3 x 10-4 
V(E) cm s-4 : 3°9 X10 3°9 X 104 
V(1 000) cm s~!) VS ~VS 
Vihg, hp, he) VS VS 
3(E) km . 61 
8(1 000) km . 4°7 X10 ~4°7 X 108 
8(hg) km . 4°3 X 1012 4°3 X10%2 
3(hp) km . 1°8 x 1016 1°8 x 1016 
3(h,) km . 2°9 x 1088 2°9 x 1018 


610 


An examination of the above two tables shows that from the point of view of 


the transmission of geomagnetic disturbances, the atmosphere may be divided into 
two parts, separated by a transition region. 


(i) At the 1ookm level the value of wr does not exceed ~ 10° while 7 has a° 
value 3 x 107, so that for waves with periods in the range listed transmission is 
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always by mode (b) of Section 3. In this mode the gas as a whole does not move 
appreciably but the magnetic field diffuses through the gas. A disturbance suffers 
a time delay in passage and, if oscillatory in form, its amplitude is reduced. 

This lower (L) or diffusion region extends from the ground up to levels of a 
few hundred kilometres depending on the value of w. The transition region 
occurs where wr ~ 7. For waves of period 1s this is between 200 and 300km; 
for waves of period 104s it is between 500 and 1ro00okm. The L region medium 
is dispersive. 

(ii) Above the transition region we have wr > 1+7 and transmission is by 
mode (c) of Section 3, that is by hydromagnetic waves travelling in the plasma 
alone, the latter moving freely through the neutral atom gas. This upper (U) 
region extends into interplanetary space. Absorption effects in this region are 
generally negligible. 

Although absorption in the U region is generally unimportant, transmission 
may be very complex. In the first place, even in a uniform medium, the O hydro- 
magnetic wave* is transmitted anisotropically. The Poynting flux, averaged over 
a cycle, flows along the magnetic field lines. Thus a disturbance having a con- 
figuration corresponding to the O wave and originating several Earth radii away 
would tend to be conducted exclusively into regions of high latitude. This may 
explain the greater amplitudes of some disturbances in polar regions and is dis- 
cussed in connection with micropulsations in Section 9. 

On the other hand the E wave travels freely in any direction. The collision 
of a solar gas cloud with those lines of force lying near the noon meridian will 
initiate an E wave which travels to all parts of the Earth, increasing the horizontal 
component of the field as observed. The delay in reaching the dark side of the 
U region corresponds to a velocity of propagation 107 — 108cms~! and a distance 
~ 2x 104km giving a time of 20— 200s. 

In addition to the peculiarities of the O and E wave in a uniform medium, one 
must also take into account the fact that the values of p, p, and Hp all vary, often 
changing substantially within a fraction of a wavelength. There should therefore 
be refraction and reflection effects including standing waves. An examination of 
these phenomena may be made by the methods of Ferraro (1954) and Ferraro & 
Plumpton (1958). 

These conclusions suggest that the analyses of Ashour & Price (1948), Sugiura 
(1950) and others might require considerable extension before they were applic- 
able to disturbances originating above a few hundred kilometres. 


6. The movements and locations of the geomagnetic current systems 


Before proceeding further it is useful to prepare a simplified model of the L 
and U atmospheric regions in the equatorial plane. This is shown schematically 
in Figure 1a, the magnetic lines of force being shown by the dashed lines. The 
limits of the U region are denoted by the planes A and B which may move but 
must at all times contain fixed quantities of gas and magnetic flux. The limits of 
the L region gas are C and G, which may not move at all under the influence of 
disturbance fields. However, the field normally contained between C and G may 
be increased or decreased by a downward or upward movement of the U region, 
which may penetrate into, or separate from, the L region. 


* The O wave has a perturbation magnetic field (H—He) perpendicular to the plane containing 
the wave normal and the steady field Ho; the E wave has its field in this plane. 
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Suppose now that equilibrium conditions are disturbed by the application of 
a force F on the upper boundary of the U region as in Figure 1b. This might 
be due to the arrival of a solar gas cloud as in the Chapman & Ferraro (1931) 
theory of the first phase of a sudden commencement (SC) storm. An electric 


current flows in the region A corresponding to the initial compression of the field 
and given quantitatively by 


curl H = 47j. (19) 


This current does not directly cause a geomagnetic disturbance because of the 
effectiveness of the U region as a shield. However, it does result in a force of 
density jAH on the gas and so causes a downward travelling hydromagnetic E 
wave which has its own current system. When the wave reaches level B it causes 
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Fic. 1.—The upper and lower atmsopheric regions. The Earth’s 
magnetic field (shown as dashed lines) is frozen into the U region but 
moves through the L region. 


the gas and field of the U region to interpenetrate the gas of L region. The field 
of L region is compressed and so the disturbance field diffuses through to the 
ground level G. 

It is the current system in the L region which is directly responsible for the 
observed geomagnetic disturbances. The Earth is effectively shielded from any 
current systems located in the U region. This result holds equally for regions 
away from the geomagnetic equator; the geometry of the field is more complicated 
but the general conclusion is unchanged. 

It has often been suggested “‘that important sources of the field of geomagnetic 
storms were located in the atmosphere” (Forbush & Vestine 1955; by atmosphere 
they mean the region below ~ 1o0okm). The present analysis shows that the 
currents directly responsible for all geomagnetic disturbances must flow below 
~ 1000km. This result might provide an argument against the atmospheric 
dynamo theory of magnetic disturbances developed by Wulf and Vestine (see 
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Vestine 1953). Observational evidence that some currents flowed in the ionosphere 
favoured a purely ionospheric origin of these currents. It would seem, however, 
that this result is explained equally well wherever the original source of the 
disturbance. 

The above simple picture of a current system propagating through the U 
region by hydromagnetic action is not, alone, sufficient to account for all observa- 
tional data. For example, Forbush & Vestine (1955) have shown that the increase 
of the horizontal component during an SC increase is anomalously large at the 
equator during the day and larger on those days when the diurnal variation Sz is 
high. This is probably due to the high conductivity of the ionosphere at the 
times and place referred to, together with an electric field system associated with 
the high-latitude disturbance. An electric field system at levels below a few 
thousand kilometres might explain other inferred current systems and also the 
aurora, 

An electric field system might be introduced near the poles, either by an 
influx of positive and negative ions along different paths or in association with a 
wave, presumably hydromagnetic. The former method meets certain difficulties 
which are discussed in the following section. An apparent objection to the latter 
lies in the fact that a plane hydromagnetic wave in the U region has a very weak 
electric field. This is because of the very high “effective” conductivity o3 in the 
plasma alone, allowing a weak electric field E to drive the current required. How- 
ever, it is well known that the effective conductivity is o3 only when the Hall current 
o2H x E/H is prevented from flowing by charge accumulations on boundaries. 
When this current may flow in a closed circuit, then the effective conductivity 
may be reduced. Waves with this characteristic may possibly have associated 


electric fields large enough to account for aurorae and the current systems of 
polar origin. Such waves comprise a twisting motion of a bundle of magnetic 
lines of force and would seem to warrant further study as a possible factor in 
magnetic disturbance theory. 


7. The penetration of the solar cloud into the Earth’s magnetic field and 
the location of the auroral zones 


In the previous sections the transmission of geomagnetic disturbances from 
several Earth radii has been briefly considered. It is likely that most of these 
disturbances originate in the region of collision between the solar gas cloud and 
the Earth’s magnetic field and a study of this collision region is an essential part of 
the transmission problem. In particular, we are concerned with the degree to 
which the solar gas cloud penetrates the Earth’s field in addition to distorting it. 
As will be evident later, the question of penetration is the central problem in the 
theory of the main phase of a storm and is highly significant in auroral theory. In 
order to avoid confusing ‘“‘penetration” with “distortion” the former is defined 
here by the co-latitude @ of the point on the Earth’s surface and in the meridian 
plane passing through the Sun, cut by the magnetic line of force which just touches 
the surface of the cloud. Lines of force lying in the same plane and cutting the 
Earth at points with smaller co-latitudes will all penetrate the cloud; those with 
larger co-latitudes will not. 

The problem of the approach of a cloud with a plane surface towards a dipole 
has been considered by Chapman & Ferraro (1931). Near the Sun-Earth line 
the dipole field is compressed and the forward motion of the gas is arrested. 
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North and south of this line the opposing force falls to zero before rising again so 
that “horns” of gas develop on the surface of the cloud, penetrating freely into 
regions of high latitude (see Chapman & Ferraro’s Figure 4). 

When this theory was developed it was not known that the Earth’s field was 
permeated by a highly conducting plasma. The presence of the plasma necessi- 
tates important changes; in particular it eliminates the horns of gas. When a 
rigid, perfectly conducting sheet approaches a dipole in free space the magnetic 
field ahead of the sheet is found by adding the fields of the original dipole and an 
image dipole behind the sheet. Components of the two fields pointing in oppo- 
site directions cancel one another. This means that at the points on the sheet 
where the original field is perpendicular to the sheet the field becomes zero and 
if the sheet were replaced by the surface of a gas cloud, horns should develop 
around these points. This means, in effect, that as a particular line of force of 
the original field becomes perpendicular to the sheet it is “‘cut’’ and instead of 
remaining a loop its ends go to infinity in the North and South directions as sketched 
by Chapman & Ferraro. 

In the real problem the Earth’s field is “frozen into” a plasma and the field 
due to currents in the cloud is “frozen out of” that plasma. If the plasma were 
perfectly conducting the separation would be complete and none of the lines of 
force of the Earth’s field would be cut. Even those lines leaving the Earth at the 
magnetic poles would remain complete loops, suffering distortion but excluding 
the solar gas. The presence of a magnetic field in the cloud before it nears the 
Earth would not affect this conclusion. Assume that initially the face of the cloud 
is a plane; ahead of the cloud a hydromagnetic compression wave is formed and 
traverses the U region as the cloud is brought to rest. In the quasi-stable situa- 
tion ensuing, the whole of U region is compressed but is not invaded by any 
external field or gas. 

Near the Sun—Earth line the forward motion ceases when the hydromagnetic 
pressure difference across the face of the cloud equals the momentum destroyed 
per second. If we can neglect gas pressure gradients then 


H2 = 8mpu2 (20) 


where H, is the strength of the distorted field just ahead of the cloud and p and u 
are the mass density and velocity of the solar plasma. This result follows whether 
or not the dipole field is frozen into a plasma. However, the physical situations 
are different in the two cases: in the former case the plasma and field are com- 
pressed while in the latter case the field is increased by the invasion of a field of 
external origin which is parallel with the dipole field. 

Away from the Sun-Earth line the different physical situations are really 
important because the oppositely directed fields, which would cancel one another 
in a vacuum, are kept apart by the plasma. Thus, cutting of the lines of force and 
penetration by the solar gas is prevented. 

The forward motion of the neutral atoms in the solar cloud is not arrested by 
the Earth’s field. They continue on into the atmosphere where some are ionized 
(they will not lose much energy in this process) and trapped in the field. These 
protons have energy up to a few times 1o*eV. 

Since the solar plasma alone, neglecting neutral atoms, has finite conduc- 
tivity (os), the compressed magnetic field will diffuse through its surface. The 
situation, neglecting distortion of the surface of the cloud away from the Sun— 
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Earth line, is shown in Figure 2. The surface is at rest in the plane x = xo with a 
particular magnetic line of force Lg touching its surface and cutting the Earth’s 
surface in co-latitude @. The direction of the field near x = xo is parallel to Oz. 
The equation giving the rate of diffusion of the field into the cloud is equation (16) 


above. If diffusion starts at time ¢ = o, then the field in the plasma is given, after 
integration, by 


H(x) = H,.exp(—703(x— xo)?/t). (21) 


The total magnetic flux through the cloud and lying between parallel planes 
y = const., separated by 1 cm, at time ¢ is 


= f H(x)dx. 
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Fic. 2.—Distortion of the lines of force of the Earth’s field by a solar 
gas cloud. 


Using equation (21) and integrating 


t \+ 
H.(—) ; (22) 
403 
Through the geometry of the undistorted dipole field of the Earth, the flux % 


may now be related to the co-latitude @ of the line of force Lg. The first relevant 
equation is 


~=M | xo 8dx (23) 


where M is the Earth’s dipole moment and ¢ the flux per cm thick sheet lying 
beyond a distance xg which, in turn, is related to the co-latitude by 


Xe = asin~26 (24) 


where a is the Earth’s radius. Combining these equations and putting 
M = 8-2x 10% and a = 6°37 x 108, we have 


sin8@ = 2-5 x 10°17H2to3". (25) 
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The value of H, is found from equation (20); for a stream of (large) proton 
density 1 00ocm~3 and velocity 108cms~! it is 0-02 gauss. The value of og is at 
least 10-8e.m.u. and if 6 = 23°, corresponding to the auroral zone, the time 
taken to penetrate is 5-5 x 108s or 18 years. The time available for penetration is 
a few hours so that according to the above analysis the flux which may penetrate 
the cloud is much less than that between the pole and the auroral zone; it corre- 
sponds to that down to co-latitude of about 6°. 

This result raises some doubt concerning Martyn’s (1951) theory of the loca- 
tion of the auroral zone. Martyn’s criterion of penetration was that the stream 
advanced into the Earth’s field until the pressure exerted by the lines of magnetic 
force equalled the pressure of the stream. In the above analysis this is a criterion 
of distortion rather than penetration of the field. 

Although the above discussion applies quantitatively only near the Sun—Earth 
line it seems applicable to theories of the direct entry of gas into auroral regions. 
Whatever path is followed by this gas it must cross a total magnetic flux given by 
equations (23) and (24) and when @ corresponds to the auroral zones (~ 23°) 
this amount is about 104 times larger than is allowed by simple diffusion theory. 

There are a number of other factors which may influence the penetration of 
the solar cloud; these have been neglected in the above discussion and in the 
earlier treatments referred to (Chapman & Ferraro 1931, Chapman 1952, Ferraro 
1952, Martyn 1951). First there are the effects of electric space-charge fields 
which will build up on the east-west boundaries of the Earth’s magnetic field 
and also perhaps in small-scale distributions due to the inherent instability of a 
medium comprising two or more separate electron streams. There is also the 
possibility of hydromagnetic turbulence near the boundary between the fast- 
moving cloud and the U region. Beyond recalling that large-scale space-charge 


fields tend to be short circuited by the highly conducting surrounding gas (Pid- 
dington 1958), these questions will not be pursued further here. The results of 
the next section suggest that some mechanism allows substantial penetration, 
but only after a period of about half a day. 


8. The main phase of an SC storm 


The main phase of an SC storm is a worldwide decrease in magnetic intensity 
lasting for a day or two. If the space outside the Earth’s surface were a vacuum 
then the decrease would be similar to that caused by an electric current flowing 
near the geomagnetic equator in a ring of radius a few Earth radii. Chapman & 
Ferraro and later Martyn (1951) have invoked such a ring current to explain 
the main phase and other storm phenomena. 

However, the space outside the Earth is not a vacuum and the presence of 
ionized gas greatly modifies the problem, as pointed out by Parker (1956). Parker 
and also Singer (1956) criticized the ring current theory on the grounds that the 
conductivity of the interplanetary gas is so high that the field due to the ring 
current would take months or years to penetrate to the Earth. Strictly speaking 
this objection is invalid because the effective electrical conductivity of the gas 
concerned is very low (Appendix and Table 1). However, because of the large 
collision period 7, the fully ionized plasma must, in this particular problem, be 
considered to move independently of the neutral atoms. The plasma alone has a 


high conductivity and is impervious to magnetic fields in the period of time 
available. 
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As seen above, the presence of this plasma necessitates changes in the theory 
of the SC; corresponding changes are required in the theory of the main phase. 
The immediate effect of a ring current is to reduce the original field just inside 
the ring and increase the field just outside. This localized effect is transmitted 
downwards to the Earth hydromagnetically, the whole of U region being lifted 
and the field in Z region then reduced. The situation is illustrated in Figure 1(c) 
where the U region has been lifted to leave a space BC between the U and L 
regions. Such a space is only formed in the simplified model; in nature there is 
no surface of discontinuous change in properties of the medium; one merges into 
the other. 

The fundamental physical feature of the main phase is an outward radial 
mechanical force which raises the U region. Such a force is not provided in 
theories where the ions and electrons spiral around the Earth in opposite directions 
(Singer 1956). The ions and electrons must move in the same direction with 
nearly the same velocity so that for a given current the centrifugal force is sufficient 
to balance the Lorentz force. This condition is satisfied in Martyn’s (1951) 
discussion of the ring current theory. There appears to be a difficulty, however, 
caused by the presence of an electrostatic field which must exist as a necessary 
part of the mechanism (to compensate for the different centrifugal forces on ions 
and electrons). This field is experienced by the gas occupying the ring volume 
prior to the entry of the solar ions whose movement provides the ring current. 
The effect of the crossed electric and magnetic fields, E and H, on the plasma is 
to accelerate it to velocity E/H; this velocity turns out to be the same as that of 
the circling solar ions (Martyn 1951). The main phase of a storm appears to 
build up gradually from zero so that initially E = o and the velocity of drift is 
zero. This would mean that individual solar ions entering the ring must share 
their kinetic energy with all other ions in the ring to provide a gradual increase 
in E and a corresponding acceleration of the gas. It is difficult to see how this 
could occur, as it would mean that the ions were giving up their energy to a steady 
magnetic field in the absence of an electric field. 

Any explanation of the main phase must depend primarily on penetration of the 
solar gas into the geomagnetic field. Only after penetration can the solar gas exert 
an outward pull on the field. There appears to be a simpler explanation of this 
pull than the delicately balanced ring current: the magnetic lines of force which 
diffuse into the solar gas may be carried away by the gas flowing past the Earth to 
provide a field something like that of Figure 3. Arguments in favour of the 
formation of such a field are similar in some ways to those given for a radial solar 
field (Piddington 1958). Once a magnetic tube of force enters the solar gas which 
is sweeping past the Earth it is frozen into that gas and drawn out beyond the 
Earth’s orbit. The additional magnetic energy so created is provided at the 
expense of the kinetic energy of the cloud. A trailing field of strength H, exerts 
a pull of H,?/8m per cm? on the receding cloud. 

On the above hypothesis the sequence of events during an SC storm (ignoring 
special polar effects) is as follows. The solar cloud envelops and sweeps past the 
Earth’s field compressing it on the solar side. Near the geomagnetic equator the 
field is everywhere nearly parallel with the surface of the cloud as in Figure 2 
and the force F = H,2/8m is transferred downwards to the surface of the Earth 
as in a purely hydrostatic problem. At the surface of the Earth it increases the 
magnetic pressure by approximately the same amount at all longitudes. As time 
passes, the solar gas penetrates the field so that its ordered motion of 10®cms~1 
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is converted, at least in part, to disordered motion. The corresponding tempera- 
ture of the protons is 4 x 107°K, so that a small amount of gas can expand the 
field until its outermost lines of force can be blown away. The lines of force so 
drawn out will cut the Earth at all longitudes because a reduction of magnetic 
pressure within a particular range of longitudes will cause an influx of field from 
outside that segment. The result is a decrease in magnetic pressure, and field 
strength, at all longitudes which, when it grows to exceed the increase due to the 
force F, constitutes the main phase of the storm. 





Fic. 3.—Distortion of the lines of force of the Earth’s field during the 
main phase of an SC storm. 


Clearly the penetration is the central problem and its mechanism is not evi- 
dent. Nevertheless it is possible to estimate, very approximately, the degree of 
penetration. Let H and h be the undisturbed field and the decrease during 
the main phase. The decrease in magnetic pressure is then everywhere 
in{ H2—(H—h)?}. At the top of the U region where the field is entirely torn 
away we have h = H = Hg say, the pressure change being H,?/87. If this is 
equated to the change at ground level, where h < H, then we have approximately 


H? = 2hH = o-66h. 


For a storm of rooy (hk = 10-% gauss), Hg = 0-026gauss, which corresponds to 
the undisturbed dipole field at a distance of z ~ 2-4 Earth radii. The lines of force 
cutting the geomagnetic equatorial plane at this distance return to the Earth at 
co-latitude given by sin 8 = (2°4)-#, @ = 40°. 

These results, which may be only rough approximations, suggest two things. 
The value of H, is much larger and the value of z much smaller than found from 
ring current theory. This is because the geomagnetic disturbance is transmitted 
through the U and L regions as a hydromagnetic pressure change rather than 
directly as an invading magnetic field of external origin. Thus the disturbance 
at a few Earth radii may be much more violent than previously believed. The 
second item of significance is the large value of 0, suggesting a mechanism of 
penetration much more efficient than that dealt with above. 

These ideas will be discussed in more detail elsewhere, together with possible 


significance in connection with Gegenschein, whistling atmospherics and cosmic 
ray variations. 
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9. Micropulsations 

If the ordinary equation for the rate of absorption of hydromagnetic waves in 
a single fluid (equation (11) above) is applied for waves of periods 100s or less 
in the E-region, they are damped out in a metre or less. This result seemed to 
provide an objection to the otherwise attractive suggestion (Dungey 1954, Lehnert 
1956, Kato & Akasofu 1956, Obayashi & Jacobs 1958 and others) that micro- 
pulsations resulted from hydromagnetic waves travelling to the Earth from inter- 
planetary space. However, according to the above theory, electromagnetic waves 
of the periods concerned do not traverse the E-region as hydromagnetic waves but 
rather diffuse through a stationary conductor without excessive losses. 

This result, if correct, would invalidate the theory of Dessler (1958) who 
assumed that the disturbances travelled as hydromagnetic waves in the ion- 
electron plasma alone. He found a minimum wave velocity at ~ 400km of 
~ 10’cms~! followed by a rapid increase in velocity at lower levels (due to the 
decreasing ion density). According to Table 2, however, the velocity is lower at 
lower levels because of the effects of neutral atoms. 

As mentioned in Section 5, some properties of micropulsations may be ex- 
plained in terms of O type hydromagnetic waves. Referring to Figure 2, the 
deformation of the Earth’s field at the plane x = x9 may propagate directly to the 
Earth as an E wave. Now consider lines of force lying in the plane x = 0. The 
gas cloud moves past the outer edges of these lines and must deform them so 
that they lie on a curved sheet whose inner portion remains in the x = o plane 
and whose outer portion is nearly parallel with the y = o plane. The main 
deformation occurs near the region where magnetic and inertial forces balance 


according to equation (20)*. The deformation field is parallel to Ox so 
that a hydromagnetic disturbance may only travel to the Earth along the magnetic 
lines of force. As Martyn (1951) has shown, reasonable values of p and u in 
equation (20) give H, corresponding to a line of force leading into the auroral 
zone. This might explain the fact that micropulsations are observed most fre- 
quently in those latitudes. 


Radiophysics Laboratory, 
C.S.I.R.O., 
Sydney, 
Australia: 
1959 January. 


APPENDIX 


The effective electrical conductivity of the outer atmosphere 
and interplanetary space 


An ionized gas permeated by a magnetic field Ho is anisotropically conducting, 
having conductivity components a9 (parallel to Ho), o; (Pedersen) and o2 (Hall). 
The “effective” conductivity, which determines the rate of diffusion of a magnetic 


* The physical concept behind this equation differs from that given by Martyn and by 
Obayashi & Jacobs (1958) who compare gas and field energy densities. It is similar to that of 
Warwick (1957), when his gas slab has thickness (h) equal to the radius of curvature (r) of the plane 
containing the line of magnetic force. 
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field through the gas, is then o3 = 01+ 22/0) (Piddington 1955, 1956; Cowling 
1956). 

In the particular case of a fully ionized gas og = oo and the effective conduc- 
tivity for hydrogen is given approximately by 


go ~ 1074732 e.m.u. (26) 


where T is the kinetic temperature. At a distance of two or three Earth radii 
where the temperature may be about 104°K the conductivity would then be 
about 10-8e.m.u. The time taken for a magnetic disturbance to penetrate to the 
Earth by a process of diffusion through a stationary medium would be some years. 

However, the gas concerned cannot be fully ionized; there must be a balance 
between the processes of ionization and recombination, with some neutral atoms 
present. It will be seen that a very small proportion of such atoms, say 1 in 1019, 
may reduce the effective conductivity below that given by equation (26). 

Neglecting the unimportant dissipative effect of electron-neutral atom collisions, 
a general expression for the effective conductivity of a gas is 


I pntH? 
PE 


; : 2 
03 90 (p+pn)? 7 


where p, pn are the mass densities of the ion plasma and neutral atom gas and 
t is a collision period given by 
3 { mimo(m, + me) \' 
T= —— a> Ps 
Spr2\ = axkT OJ 





(28) 


Here r is the collision distance for heavy ions and neutral atoms and mj, m2 the 
masses of those particles. The physical significance of the unusual collision 
period chosen is discussed below (see also Piddington 1957, Appendix I). 

Accurate values of r for the conditions concerned are not available but a 
reasonable estimate is probably given by (Cowling 1945): 


/200\* 
r= 10°*K( =) ’ (29) 


where K ~ 6 in the Earth’s true atmosphere and K ~ 1 at a few Earth radii 


where the gas is presumably mainly hydrogen. Combining equations (26) to (29) 
we have 


I 1014 8-1 x 10°3N,,H2 


og T32 N(N+N,)®AtK? (30) 
where N, N, are the number densities of ions and neutral atoms and A their 
atomic weight, assumed equal. 

The values of the various quantities determining o3 are discussed in Section 4. 
For the present purposes a set of likely values, corresponding to a distance of a 
few Earth radii are chosen as follows: N = 103cm-3, Ny = 10-3cm-3, A = K = 1; 
neglecting the small term in oo, equation (30) then reduces to 


o3 = 1'2x 107222, (31) 


The value of H in the equatorial plane is about 0-332~3 where z is the distance in 
units of Earth radii. At say 5 Earth radii we have og = 1-7 x 107!’ e.m.u.; the 
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value of H is of the order of observed geomagnetic fluctuations. The value of o3 
decreases as the Earth is approached until the ionospheric E region is reached, 
where it rises to a maximum value of about 10-!%e.m.u. (Baker & Martyn 1953). 
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Travel-times of PKP up to 145° 


B. A Bolt 
(Received 1959 March 16) 


Summary 


Examination is made of some goo PKP readings in the range 
110° < A < 145°, as recorded in the International Seismological Sum- 
mary from 24 southern earthquakes. This provides new evidence 
that the J.-B. times for the DE branch before 140° need an increase 
of about a second. There is no evidence in this range of significant 
differences between PKP times to Europe and North America. 

In the range 135° < A < 145°, a concentration of PKP observa- 
tions preceding the DE branch is found fitting the form 


t = 19min 29:18 + 2°69(A — 141°°6). 


This possibly corresponds to the BC branch of PKP and its extension 
for the diffracted waves associated with the caustic at B. 


1. Introduction 


One significant difference between the independent solutions of Jeffreys 
(1939a) and Gutenberg (1951) for the variation of P seismic velocity with depth 
in the Earth occurs near 5000km. Jeffreys’s solution has a negative P velocity 
gradient between depths of 4982km and 5121km. At 5121km there is a dis- 
continuous velocity increase marking the boundary of the inner core. Gutenberg 
however finds that no velocity solution consistent with his 1951 travel-time data 
gives ‘‘a decrease in velocity or a strictly vertical segment”’. 

At depths near 5000km both solutions depend considerably upon observed 
travel-times of PKP phases, particularly the phase PKIKP which enters the 
inner core and is of P type throughout its path. Recently, Bullen & Burke-Gaffney 
(1958) have shown that the 1954 hydrogen bomb explosions yielded some impor- 
tant additional data on PKP travel-times. The purpose of the present work is to 
present further information, obtained with satisfactory precision from PKP 
observations from natural earthquakes. 

The Jeffreys-Bullen PKIKP travel-times depend largely on PKP readings 
from the Solomon Islands earthquake of 1932 January 9 (Jeffreys 1939b). They 
were subsequently checked from observations of the Celebes Sea earthquake of 
1934 June 29 (Jeffreys 1942). The present paper presents an analysis of nearly 
goo PKP observations over the range 110° < A < 145° from 24 selected southern 
earthquakes, and provides some support for the results of Jeffreys. 


2. Theoretical PKP travel-time curve 


The form of the PKP travel-time curve is shown in Figure 1. It is divided 
into four branches by the three cusps B, C and D. 
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A ray which just grazes the outer core is refracted to A at A x 180°; steeper 
rays give in turn the branches AB and BC until at C the ray is tangent to the inner 
core. The cusp B at A x 142° is similar to an optical caustic. The receding 


branch CD corresponds to rays reflected by the inner core, and DEF to the PKIKP 
rays refracted through the inner core. 


A 


Time (t) 








——* Distance (A) 
Fic. 1.—Theoretical travel-time curve for the core phase PKP 
(Jeffreys & Bullen 1948, Figure 2). 


3- Reduction of readings 

Three criteria were used for selection of earthquakes: 

(a) occurrence in the Celebes to New Hebrides region over the period 1937 to 
1941, 

(b) at least 12 readings of PKP for A < 145° given in the International Seis- 
mological Summary ; and 


(c) no evidence of multiplicity. 
Only the 24 earthquakes listed in Table 1 were satisfactory. Six others were 


rejected as multiple after examination of the Riverview records and the residuals at 
all stations. 


Table 1 
The selected earthquakes 
1937 January 23 1938 May 30 1939 October 17 
January 25 June 9 1940 February 20 
July 2 August 31 March 21 
August 5 October 20 May 28 
September 15 1939 January 30 September 26 
September 27 February 3 1941 September 4 
1938 February 1 April 5 September 9 
May 12 August 18 November 8 
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The J.S.S. solutions were altered in six cases by the method of Jeffreys 


(1939c, p. 131) and retained in eighteen. If x and y are the north and east dis- 
placements to the J.S.S. epicentre, the altered solutions are: 


Revised Revised 
x y focal-depth time of origin 
km h mins 
1937 Sept. 15 _ _ 78+8 12 35°1+0°7 
1939 Feb. 3 —0°31 40°03 +0:06 +0:05 —_ 5 19°9+0°4 
1940 Sept. 26 +0:14+0-01 —0'41+0°19 150+3 3 33°8+0°4 
1941 Sept. 4 —_— — 88+5 10 46-0+0°6 
1941 Sept. 9 —o'58+0-02 —0°14 10°24 — | 47°240°5 
1941 Nov. 8 +0°53 10°03 +0:03+0°10 86+3 23 31°71+0°6 


Using the adopted solutions of all 24 earthquakes, the ratio of the number of P 
residuals between + 4s (inclusive) to the total number recorded at each of the 
following near stations is: Brisbane 0-82, Riverview 1-00, Melbourne 0°68, Ade- 
laide 0-70, Perth 0-43, Christchurch 0-92, Apia 0-83, Manila 0-85. These stations 
have sufficiently different azimuths for the ratios to imply that significant E-W 
or N-—S systematic errors in the adopted epicentres are unlikely. 

For Pacific sources there is evidence (Jeffreys 1954, Burke-Gaffney & Bullen 
1958, Carder & Bailey 1958) that the J—B. P tables need corrections of about 
—2s. Thus, for each shock, an adjustment equal to the mean P residual for 
80°< A < 105° was applied to the observed PKP times; this also partly removed 
discrepancies arising from incorrect depth assessments. In general, the adjust- 
ments were less than 1s, the maximum (occurring with 1941 November 8) being 
2*1+0°7S. 

PKP residuals for all stations with vertical component instruments were calcu- 
lated using the J.-B. tables for the branch DEF (adjusted for ellipticity). All 
phases within 20s of these J.-B. times were included. The summaries of the 
frequencies of the residuals for European and North American stations are shown 
in Tables 2 and 3. (For statistical convenience the length of the range interval 
has been twice reduced.) 

Tables 2 and 3 exhibit the following prominent features. 

(a) For 110° < A < 140°, there is a main concentration of observations, 
taken to refer to the DE branch of PKP, with its mode between 1s and 2s. As 
found by Miss Lehmann in her study (1958) of the Celebes shock of 1948 March 1, 
there was no concentration for A < 110°. 

(b) For 135° < A < 145°, there is a separate distribution of residuals from 
observations which clearly occur before the main concentration noted in (a). 


These are possibly associated with the branch BC of PKP and will be treated in 
Section 5. 


4. Travel-times for the DE branch (PKIKP) 


Before combining the summaries in Tables 2 and 3, a test was made to deter- 
mine whether PKIKP travel-times to the same distance along European and 
North American paths are equal. The means in Tables 2 and 3 were formed after 
allotting unit weight to each observation in the square brackets. (Detailed analysis 
indicates that this gives satisfactory estimates of the mean and its standard error 
in each range.) The most significant difference between the means is — 1-0 + 0°798 
for 132°-6 < A < 135°-0, and for the 8 differences, x? = 8-8. Allowing for 
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uncertainty in the epicentres it follows that PKIKP travel-times for the two paths 
are equal within about 1s. 

The two sets of PKIKP residuals in square brackets may now be combined. 
Following the statistical method of Jeffreys (1939c, p. 162), the law of error for 
the residuals x; is taken to be 


S(xi) = 7*(1 — m)hexp{ — h2(x4— «)?} + mg(xi), 


where h is the precision constant, « the weighted mean, and g(x;) a slowly varying 
function estimated using the method of uniform reduction. ‘Taking six ranges at 
intervals of 5°, assumed means of 1,2,1,1,1,28 respectively were subtracted and 
the frequencies then totalled to give the summary in Table 4. 


Table 4 


Observed and calculated frequencies of the combined PKIKP residuals for 
110° < A < 140°-0 
x, (sec) —-8 — —6 5s —4 
Observed 18 
Calc. I 9 20 
O0-C —2 -2 
Calc. II 15 


x, (sec) I 3 4 
Observed 68 38 26 
76 37 20 

—8 I 6 

64 33 21 


The mode of the distribution is approximately zero so an estimate of the standard 
error, c, can be obtained from the ratio of the number of residuals between + 4°58 
to the number between +2-5s. The probability that x; is between +a is erf ha 
(Jeffreys 1939¢, p. 61) so that 


erf4-5h 421 


= = 1°32, 
erf2-5h 318 
where 2h?6? = 1. 

The solution is h = 0-30 so that computing 83 exp(—A?x;2) we obtain the 
frequencies called Calc. I in Table 4. The differences (O—C) suggest that the 
asymmetry can be dealt with by a uniform reduction of 4+ 42. This gives Calc. II 
in Table 4, with h = 0-29 and o = 2-48. Adopting this approximation, the mean 
for each range, shown in Table 5, was then calculated using appropriate weights. 


Table 5 


Means and unsmoothed travel-times for PKIKP 


Mean Mean Unsmoothed 
distance residual travel-times 
(°) s min s 
112°6 +1°49+0°55 18 39°8 
117°6 +1°98+0-41 50-0 
122°6 +1°4440°25 59°1 
127°6 +1°11 +0°23 19) | 85 
132°6 +1°41+0°24 18-3 
137°6 +211 +0°29 28-2 
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These means were used as corrections to the J.-B. PKIKP tables to give the 
‘“‘Unsmoothed travel-times” in Table 5. ‘These were first fitted to the form 


t= 19min 3s+9x+a+bx+cx?, 


where 5x = A—125°. The formal solution is 
a = 0°57+0°17; b = 0°53 + 0°03; ¢ = 0°09 + 0°02. (1) 


By theory, if all included observations are PKIKP, DE is curved downwards 
and ¢ is negative, so that this solution presents some difficulties in interpretation. 
(The upward curvature may disappear if the present data were combined with 
similar observations for A > 145°.) Since the quadratic term in (1) contributes 


a maximum of only o-8s at 110° and 140°, we alternatively fit the times to a 
linear form, yielding 


t = 19min 3°78 + 1°92(A— 125°)s, (2) 
with o = o-6s. 

For 110° < A < 140°, the mean difference between (2) and the J.-B. tables for 
the DE branch of PKP is +1°5s. This agrees with the mean PKIKP residual of 


1:28 for A < 140° in Jeffreys’s preferred solution (1942, p. 36) of the 1934 June 29 
earthquake. 


5. Travel-times for 135° < A < 145° for the early observations 


Some of the residuals for A > 140° may belong to PKIKP. ‘To reduce the 
effect of these, consideration was limited to the residuals in the curly brackets in 
Tables 2 and 3. After subtraction of rough means in each range a total distribution 
was formed which, apart from a slight scarcity of observations in the central 
frequency, resembled a normal one with o = 1°9, superimposed on a uniform 
background of 2s. 

Proceeding as in Section 4, seven weighted means were used to obtain, for 
135°:0 < A < 145°-0, the linear form 


t = 19 min(2g9°1 + 0°12)s+(2-69 + 0°03 )(A— 141 °-6)s, (3) 
with o = 1°19s. 

For 143°-0 < A < 145°-o, the J.-B. times for the branch BC are within 
1s of the times to the same distance given by (3) and the gradient of (3) equals the 
J.-B. gradient near 145°. This suggests that (3) corresponds to the BC branch 
of PKP and its extension for diffracted waves. 

It is of interest that (3) agrees with Jeffreys’s preferred solution of the 1934 
June 29 shock. This gives residuals against (3) of +1, 0, and os at Pittsburgh 
(140°*8), Oak Ridge (142 °-2) and Fordham (142 °-g), respectively. 


6. Discussion 


The equations (2) and (3) are compatible with PKP observations from four 1954 
hydrogen bomb explosions (Burke-Gaffney & Bullen 1958). (‘The mean P residual 
calculated from data given by Bullen & Burke-Gaffney (1958), is 8 = —2°7+0-3s 
from 17 observations with 69° < A < 97°.) 

First, allowing for 5, four readings at Stuttgart at 116°-2 give a mean residual 
of +0°2s against equation (2). Secondly, the small first arrivals at Pretoria, 
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Kimberley and Tamanrasset (Bullen & Burke-Gaffney Joc. cit.) give the following 


mean residuals against equation (3) after the adjustment 6 is applied. (m is the 
number of observations.) 


Distance 


Pretoria 137°3 
Kimberley 139°5 4 


‘Tamanrasset 140°4 4 +1°9 


There is a suggestion of a systematic deviation but agreement between the 
observations and (3) is well within the uncertainties. (Bullen & Burke-Gaffney 
remark that the uncertainty of the Tamanrasset readings is about 3s.) 

Recent investigations of the PKP phase have also been made by Gutenberg 
(1951, 1957, 1958a, 1958b) and Denson (1952). They show distributions of early 
short-period PKP phases for A < 145° which are considered by Gutenberg to 
be a consequence of dispersion near the boundary of the inner core. The present 
data, which take no account of wave periods, provide no evidence of multiplication 
of PKP branches. The derived travel-time curves (2) and (3) have very similar 
gradients to the J.-B. DE and BC curves respectively, and the standard errors 
have magnitudes usual in this type of study (Jeffreys 1939b). The results may be 
explained in terms of refraction at the boundary of the inner core and diffraction 


near the PKP caustic at B. They support Jeffreys’s suggestion (1942) that the 
radius of the inner core derived from the J.-B. tables (~ 1250km) needs a slight 
reduction. 

I am indebted to the late Fr. T. N. Burke-Gaffney, previously Director of 
Riverview Observatory, for the loan of seismograms. 


University of Sydney, 
Australia: 
1959 March 10. 
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The External Gravity Field of a Rotating Spheroid 
to the order of e’ 


A. H. Cook 
(Received 1959 March 17 ) 


Summary 


The theory of the external gravity field of a rotating body bounded 
by a spheroidal equipotential surface is carried to terms of the order 
of &. The potential is expanded in functions that are orthogonal in 
the space outside an oblate spheroid. It is found that terms of order e* 
in the gravity formula and in the relation between J, e and m are mostly 
very small. 


1. The problems 


The theory of the variation of gravity over a rotating spheroid on which the 
International Gravity Formula is based is taken as far as terms of order e2 where 
e is the ellipticity of the surface. The relation between the constants J, e and m, 
which connects the effect of the Earth on the motion of the Moon and artificial 
satellites with the variation of gravity at the surface of the Earth, is also carried as 
far as terms of order e?. Hitherto this approximation has been quite adequate 
but it is now clear that the potential accuracy of observations on artificial satellites 
requires that terms of order e? should be examined. The same order of approxi- 
mation will be required if the distribution of gravity observations over the surface 
of the Earth should be improved so that the amplitudes of the lower order harmonic 
terms can be determined to 1 mgal. This is by no means unlikely. Accordingly, 
in this paper the formula for surface gravity and the relation between J, e and m 
are taken to order e*. The Earth is not, of course, an exact spheroid but its external 
gravity field departs only very slightly from that of an exact spheroid, and so the 
theory of the external field of an exact spheroid provides an excellent standard of 
comparison for the field of the Earth itself. 

The formulae have in the past been derived from an expansion of the potential 
outside a spheroid in solid spherical harmonics but these functions are not ortho- 
gonal for the domain exterior to a spheroid and the nature of the approximations 
involved would require examination. As is well known, the theory of surface 
gravity has also been developed (Pizetti 1894) using a class of functions, elliptic 
functions, that are orthogonal outside an ellipsoid, and the formulae, which are 
exact, may be expanded to terms in e*. However, this method is less convenient 
for the study of the external field and since the spherical harmonic method for 
finding the relation between J, e and m also requires investigation, it is convenient 
to use instead a second class of functions (Spheroidal harmonics) involving Legendre 
functions, orthogonal outside an oblate spheroid. (See Jung 1956). As a check 
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on the algebra the gravity formula is established by the use of this second class of 
functions as well as from Pizetti’s formula. An incidental advantage of using spher- 
oidal harmonics is that the reason for the existence of a fourth harmonic term 
(and higher ones) in the spherical harmonic expansion of gravity outside a spheroid 
is clarified. Mineo (1930, 1933) has discussed the accuracy of Clairaut’s and 
Stokes’s theorems when derived by the use of spherical harmonics. 

The plan of this paper is as follows: 

The coordinate system and some necessary formulae are set out in Section 2. 
Section 3 contains the derivation of the gravity formula over a spheroid from the 
sole condition that it is an equipotential surface of the sum of the gravitational and 
rotational potentials and in Section 4 the relation between the ellipticity e and the 
mass and moments of inertia of the rotating spheroid is obtained by direct inte- 
gration without recourse to MacCullagh’s theorem. The expansion of the potential 
at external points, required in the theory of the motion of satellites, is derived in 
Section 5. 

It is hoped in a later paper to derive the normal gradient of gravity and to 
extend Stokes’s theorem for the determination of the elevation of the geoid from 
gravity anomalies to terms of order e. 


2. The coordinate system 
We take orthogonal coordinates (n, 9, A) defined in terms of Cartesian co- 
ordinates as follows: 


x = ccoshy sin 9 cos A, 
y = ccoshy sin $ sind, 
z = csinhy cos 9. 


We assume that there is rotational symmetry about the axis of rotation 2. 

The surfaces » = constant are confocal spheroids and the surfaces $ = con- 
stant are confocal hyperboloids of one sheet. The surfaces A = constant are 
meridional planes. is the eccentric angle and A the longitude. 

We suppose that all matter is contained within the surface 7 = yo and that 
this surface is an equipotential surface for the geopotential (surfaces with other 
values of 7 are not equipotentials). 


The major and minor axes of the spheroid bounded by this surface are respec- 
tively: 


c cosh no 
and 


b = csinh yp. 
If e = (a—b)/a, it is readily shown that 


cosh yo = {e(2—e)}+ 
sinh yo = (1—e){e(2—e)}-+. 


and 


For this coordinate system, the line element ds is given by: 


ds? = hy? dy? + ho? d3* + hg? dd, 
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hy? = c*(sinh?y + cos*$), 


ho? = c*(sinh?» + cos?) 
and 


hs? = c® cosh?» sin? $. 
The surface element is 
dS = hohgd% dr 
= ¢c*(sinh*» + cos?$)* cosh n sin $ d9 dA. 
The element of volume is 
dr = hyhoh3 dnd da 
= ¢(sinh?y + cos?) cosh n sin $ dn d9 da. 
The gradient of a function which is constant on the surface 7 is 0/h,@n, that is 
I 7 
c(sinh?» + cos?$)# én 
The expansion of (sinh2y9+cos?%)-* will often be required: 





2 4 6 
(sinh®y + cos%d)-* = I , _1 cos > ay cos*} a vo $ 
sinh no 2 sinh?y 8 sinh*yo 16 sinh®y 


to order e°, 





that is 
t 3 P 3 5 
- {1 —e(1 +—e + 2e?)cos?} + -e2(1 + 3) cost ——e3 cos*H}. (2.2) 
sinh 7 2 2 2 


Finally, we relate $ to the geographical latitude ¢. In the meridian A = 0 


Ox jdz 
cotd = — / — evaluated at » = no 
Onl On 


= (1—e) tan 9. 


Then 
cos?} = {1 —(2e—e?) cos*¢ — 4e? sin?¢ cos*} sin*¢. 


In terms of geographical latitude, (sinh? +cos?9)-* is equal to 


- {1 —e sin®¢ + fe? sin22¢ + pe? sin*¢ sin?2¢}. 
sinh no 


3. The gravity formula 


The orthogonal solutions of Laplace’s equation for the space external to the 
surface » = no which are symmetrical about the z-axis are 


Q,(i sinh n) 
~** O,(é sinh no) 


* The minus sign preserves consistency with other expressions for the gravitational potential. 


P,(cos 9),* 
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(Hobson 1931, p. 431) where the Legendre functions of the second kind are defined 
by 


+1 
_ 1 (Pale) 
Onn) = = f- du. 


—u 
-1 

On the surface yo the gravitational potential, V, of the mass inside no is accord- 
ingly 


— Xvn Pr(cos 9), 


while the rotational potential is 
4w*c? cosh2no (1 — cos? $); 


w is the angular velocity of the Earth about the z-axis. 
Since the sum of the gravitational and rotational potentials has to be constant 
over the surface 7 = 7, it follows that 


v2 = 4w*c? cosh2y9 = 4w* a? (3-1) 


while all the other coefficients except vo vanish. 

It may be noted here that when expressed in this coordinate system the 
gravitational potential over a spheroidal equipotential surface has no fourth 
harmonic term although it does have one when expressed in spherical polar co- 
ordinates. The fourth harmonic in the potential of a spheroid is thus in a sense 
fictitious since it may be removed by a suitable coordinate transformation. 

The total potential U on or outside the surface no may then be written as 

Oot sinhyn) wa? QOo(i sinh») 


“Saas ma 


The value of gravity, g, is the modulus of the derivative of U normal to the 
surface yo, that is 


i cosh oU ae 
c(sinh?y + cos?) (7 sinh ») +3 





evaluated at » = 7. 
It is convenient first to determine the constant vo. By Gauss’s theorem 


| era V .dS = 4nf M, 
Ss 


where M is the mass of the spheroid, f is the Newtonian constant and the integral 
is taken over the surface yo. 
Now 


hohg d9 dr eV 
c(sinh?y + cos?d)* dn 


= ci cosh*y sin $d9 dA 


grad V .dS = 





a(i sinh n) 
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The integral with respect to A gives a factor 27 while that with respect to 9 
eliminates the term with P2(cos 9). 


Hence 


| smav. dS = weitediat. : 
0 Ine 
s 


a dQo(i sinh ») 
Go =~ aisinh) 


ar n! I N+1 n+2 2n+3 
ak sececilnadaetinptancipehines stideias . — 
Qn) 2™*lin+4)! pm F( : + 


2.2 2 
where F is the hypergeometric function (Hobson, 1931, pp. 14 and 53). 
Hence 


I I 1 I 
O0~ o(8+ + at pa) 
and 
Qo’ = 1/cosh?. 
Accordingly 


{ grad V .dS = 4nvpic QOo(i sinh yo) 
Ss 


pf, Sea 
v = — ——Qoli sinh yo) 
ic 


{M I ig I I ) 
=——_I I — - 
b 3sinh?y%  s5sinh*n9 7 sinh®y 





f{M i.2 2 
(1 +ie+ 42), 
a 15 
The product . ™ 


Q'o 
“o- 


Qo 7, 
will be required; it is 
M 1 
fF nth G5) 
a isinhno 
The ratio Q2'/Qz2 is also needed; it may be obtained from the hypergeometric 
form of Q2 and is 
3 


8 20 
— (: —-e+ od to order e?. (3.6) 
i sinh no 7 147 
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The value of gravity at » = no is then the modulus of 








i cosh no [= I—e wa? 8 20 ) | 
x —-——-__————_ _ 


. xi — — See Sy ae 
c(sinh?o + cos?>)* a isinhyg sinh \ 7 t 2 


I 
-— — wc? cosh no sinh no x (1 — cos?$), 


c(sinh2no + cos?$)t 
that is 





I ] 5 17 
- - ~w —-e+ —#) + 
¢(sinh?n + cos?9)+ — a 
2a2/ 
+ x (5 — ~¢ — =e) cos*9 |. 
7 


2 49 
on (3-7) 
Ihe value of g at the equator, where cos $ = 0, is 


1 7a 3 


= —_—_ [= - 3 eta ( _3, - -¢)| (3.8) 
csinhyol a 2 \ 7 44! 


] 2 2 
= — — Sate (: +-e+ 52). (3-8a) 
ab 2 7 441 ) 


Now let w?a9(1—e)/fM be denoted, as is usual, by m and change from $ to ¢. 
With these substitutions, 


H. — + 3em + =m?) sin) — 


e I 
+—+ mn ) sin?2p — > etm sin*¢ sin22¢d | x 
2 2 


14 


, -. i a 
x (2 —esin?d+ ynat + =n sin®24). 


Thus, finally we have for the value of gravity over the spheroid, the following 
expression carried to the third order in e: 





| 5 1 
g.[1+(2m— e+ Smt — Lem -- = etm — Sem? + * mt) sind + 
2 4 14 98 4 8 


I 17 15 I 5 
+ gt (e—sm+ —em — — mt) sin22¢ — ;¢(= m— e) sin?d sintap |. 
8 7 2 8 \2 


} 


(3-9) | 





To the second order in e and m this agrees with the usual gravity formula 
(Jeffreys 1952, p. 131). 
At the poles, where cos $ = 1, the value of gravity is 


w= oole-9) 





External gravity field of a rotating spheroid 
to second order. Hence 
Sp & 5 


oP all 
b a 2 


(Clairaut’s theorem). 
The following formula will be found to be exact: 
2 
= + __ 4nfp— 20" (3-11) 
a b 


where f is the mean density of the matter in the spheroid (Pizetti’s theorem). 

It may be noted that although the potential of a spheroid contains no harmonics 
higher than the second when expressed in (n, 9, A) coordinates, the gravity formula 
always contains such terms because of the factor (sinh? + cos?9)-*. 


Comparison with Pizetti’s theory 


Pizetti’s result is expressed in the following formula 


1 +(8—e—eB) sin? 


Be ft —e(2—e) sin2d]* 


where 


M 
ge = — tetas + C10) 


and £, which is (gp—ge)/ge, is given by 
2w*a(1 —e) Cle) 


£e 


Bte= 


C(e) is the function 


18 125 . 
(1 +—e+ —e), 
35 245 


Consider first the value of ge. By equation (3.8a) this is 


which is identical with Pizetti’s result. 
By Pizetti’s formula, 


{1 —e(2—e) sin} + 
= 14+¢(1+e+e?) sind — $e°(3 + 2e) sin?2¢ — ge8 sin®¢ sin?2¢ 





1+B sin*d 
&e 


+ ~e(e — 5m+ tT em _ *S m2) sin? 22¢ 
8 7 2 


~ a = m— e) sin?¢ sin?2p 
8 \2 


in agreement with equation (3.9) 


4. The relation between e, m and ] 


MacCullagh’s theorem shows that in spherical polar coordinates the external 
potential of a body symmetrical about its polar axis may be written in the form 
2Ja* 


Ta tot “ef wlcns 8) + .. 
J 378 


where the coefficient J is equal to (C — A)/(Ma?), C being the moment of inertia 
about the polar axis and A that about any axis in the equatorial plane. 

This result might be used to relate e, m and J by transforming to spheroidal 
coordinates and identifying the coefficient v2 of the preceding formulae. However, 
the convergence of this particular form of MacCullagh’s theorem on the boundary 
of the body requires discussion since the conditions for an expansion in spherical 
harmonics are only fulfilled on the surface of a sphere. 

The relation between J, e and m will therefore be obtained directly without 
recourse to MacCullagh’s theorem. 

If d is the distance between two points with spheroidal coordinates (n’, 9’, A’) 
and (no, 4, A) where yo > »’, then (Hobson 1931, p. 430) 


—-=1 P (2m + 1)| Pa(cos 9)Px(cos 9)Qn(t sinh no)Pa(z sinh y’) + 2 S (-—)™x 
n=l 
7 - (n—m)! 


et ) P,™(cos 9)Pn™(cos ¥’)On™(i sinh no)Pn™(i sinh y' cos m(A— x)| : 


This formula can be applied over the surface of a spheroid where the expansion 


in solid spherical harmonics does not necessarily converge. 
The potential of the body at a point on 70 is 


fF 


and if the body has rotational symmetry, as we assume throughout, then this is 
equal to 


fi > (an+1)Pa(cos 9)Qn(é sinh yo) | dr pPa(cos 9’)Pq(i sinh 7’). 
. n 
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The coefficients of the successive terms in the expansion of V on the surface 
no are thus 


L (2n+ 1)Qp(i sinh m9) | dr pPa(cos 9’)P,(i sinh 7’). 


For n = 0, the coefficient is 


fer os 
—— .1Qo(i sinh no) 
c 
as found previously. 
For n = 2, the coefficient ve is 


sf 'O2(i sinh no) { dr p Po(cos 9’) P2{i sinh’). 


oy 
c 


To interpret this we require an expression for C— A. Taking A to be the mo- 
ment about the x-axis, 


C= [ aro(st+98) 
As | aro(y2+a 


dr = hyhohg dn’ d¥' dd’, 
x24? = c? cosh*y’ sin? 9’ 
and 
y?+22 = c%(coshn’ sin?9’ sin?A’ + sinh?n’ cos? 9’), 
We write hyheh3 = h® and then, carrying out the integration with respect to 2’, 
we find that 
C-A=cCan | dn | d9' hp(cosh®n sin? 9’ — 2 sinh®n’ cos? 9’). 


But cosh?y’ sin?9’ — 2 sinh?n’ cos?’ is equal to 


#{1 + 2P2(i sinh n’)Po(cos 9’)} 


» | dn! | a9’ Wp. 


eM 2 v2 
C-A = — 


3 15° ifQe(sinhy) 
J is defined as }(C—A)/(Ma*). Then since v2 = 4w*a?, 


2 I 2 wa2 
- J] Ma = -Me + —&—_______., 
3 3 15 3¢f Qo sinh yo) 


and M is equal to 


Hence 








i Q2(i sinh m9) = 


and accordingly 


sinh no )( 2.) -] 
= ———__—_— ——— }{1— i. + " 
2 cosh? re —e) al coshno 


I 


that is 





eh a | 2 II 
e(1 _ ~e] _ s(x —-e+ —¢) to order e°. 
- a 7 49 





This result agrees with that hitherto established to order e?. 


5. The gravitational potential as a function of spherical polar coordinates 


For use in the theory of the orbits of satellites, the gravitational potential is 
required as a function of spherical polar coordinates. The transformation from 
spheroidal coordinates (yn, 9) to spherical polars (r, @) is first required. 

We have 


r? = c® (sinh®y + sin?) 
and 
cos 8 = (c/r) sinh y cos 9. 
Put sinh » = s~!, o = c/r and p = a/r so that o? = e(2—e)p?. 
Then 
s4 cos? 6+ s*(1 — 0?)— 0? = 0, 
whence 


52 


— = 1+0°(1—cos"@) + o(1 — 3 cos?@+ 2 cos4#) + 
ot 


+ 08(1 — 6 cos*@ + 10 cos*@— 5 cos®@), 
The following are also required: 


= 1+407(1—cos?@) + $04(3 — 10 cos?6+ 7 cos4#) + 
+ #05(5 — 35 cos? + 63 cos*# — 33 cos); 


1 + $0°(1 — cos®6) + go4(5 — 14 cos?@ +9 cos*#). 
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32 
cos? = — cos?6. 
o2 


Let the gravitational potential V be written as —(Vo+ V2) where 


Qofi sinh ) 
ie = 
Qo(i sinh no) 
Q2(i sinh ») 
© Rs. dower. & 
Q2(i sinh yo) 
We take Vo first. From equation 3.4, 


Ve = 42 P2(cos 9). 


{mM 
vw = ——— Qoli sinh no) 
ic 


pee 
= —~—Qo(i sinh») 
ic 


{M I . I I 
= — {-——— ———— 
c sinh» 3sinh*y ssinh*y 7 sinh® 


s s2 s4 s6 
= fm=|1-— += - =]. 
c 3 5 7 


With the expansions just found 


= “fs + 407(1 — cos*6) + 404(3 — 10 cos"6+7 cos4#) + 4,085 — 35 cos?é+ 
r 
+ 63 cos4@— 33 cos®)} 
and 
s? 1 I 
I—-—+—— — = 1 — -o? — — o(2—5 cos*6)— 
S: Sc 3 15 


I 
— —— 0(8 — 63 cos*@ +70 cos‘@). 
105 


pete I- ~o2P2(cos 6)+ ~o4P,(cos 6) - “6 P¢(cos 6). 
f{M 3 5 7 


To obtain V2 we require the ratio 


Q2(i sinh ») 


i sinh 
which is Ort "») 
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where « = isinh y, wo = i sinh 70, 


f 3 
(@) = -G-er Gore 
be 
and (s/c)® has been given above. 


2 
F(3,2;2;uo-*) = pnt agin e 
ss 7 7 


' 6 
F(3,a;25y-2 = ~~ aoe 
- 2 7 7 6 ; 
= 1 — -o? — —o(1—6 cos? 9). 
7 7 
Combining these results, ; 


h 
(2) Si Qoe(i sinh») a -(1 204 Net) x 
1—e Qo(i sinh yo) 
x I, + a i cos? @} + — Y oM(as — 174 cos? 6+ 189 cos*é@ )). 
\ ee, 56 
Using the expression for cos $ given above 


P2(cos $) = $(3 cos?6— 1) + $0%(1 — cos?@) cos?@ 


+ §o4(1 — 3 cos*@ + 2 cos*#) cos?@. 
Hence, 
y. I fM — m Qe» o(z sinh n) 
fa t= 


——_—_—— P,(cos 9) 
3a I1-e e Ox i ee 


6 
x | Pa(cos 0) — —o* P,(cos @) + > GAPa(cos 6))}. 
7 7 
The complete expression for the potential is thus 


els 


ot ( ) ‘m(1 ~ = + — cot) fPekeos 8) + 


+ o? ial _ =(“) ‘m(1 = —- -)| x P4(cos 6) — al = — (2) m m\ Pa(cos | 
5 7 \7 
Since o2 = e(2—e)(a/r)?, we bis finally, 
M 
vy wt! ~| 


where 


eee =(; ). J P2(cos 6) + —(2 *)' D P4(cos 8) + | : )"JoPa(cos a). 


I ei). 
~ 1-=e+ a), 
7 


(5.1a) 
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as found in Section 4; 


oslo -Se)-m(e—39)] 


which again agrees to order e? with Darwin’s result; and 


Je = ——-%(6e— sm). (5.10) 
21 


The usual notation in which — 3 J is the coefficient of the second harmonic and 
3sD is the coefficient of the fourth becomes inconvenient when other harmonics 


have to be considered and it is therefore proposed that the potential should be 
written as 


f(s - Daa : )" Pa(cos.)}. 


In this notation, 





nd ee ele | 
tC ae ee 


= e?(6e — 5m). 
21 





6. Formulae in terms of mean radius and gravity 


de Sitter pointed out some time ago the advantages of using the mean radius 
and gravity in place of the equatorial radius and gravity on the equator as parameters 
in the expressions for the external potential and surface gravity. The point is 
that most measurements of geodetic arcs and of gravity are made in middle lati- 
tudes so that it is the mean values that are most directly determined by observation 
and not the equatorial values. The force of this argument is greater nowadays 
when the distance of the Moon is being observed by radar from stations in middle 
latitudes. Accordingly, the principal formulae are re-cast in terms of these para- 
meters. 

For this purpose, the exact specification of the mean must be considered. The 
radius and gravity may be expressed as functions of eccentric angle 9, of geo- 
centric co-latitude 6, or of geographical latitude ¢, and the mean values of these 
expressions differ in the terms of order e? and higher. For practical purposes it 
may be convenient to take means over geographical latitude but for theoretical 
purposes it seems preferable to take means over the angle 9; this is done here. 

The radius vector r of a point on the surface of the spheroid 7 = no is given by 

r2 . 
a cosh? — cos? 3 
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so that, remembering that ¢ cosh no = a, 


cos? cos’ $ 
2 cosh2y9 8 cosh4ng 16 cosh®y 


On expressing this in terms of P2(cos $), P4(cos $) and P¢(cos 9) and using the 
property that the mean values of Legendre functions are zero, it follows that 
r I I I 
—-at—--e+ —0+—A, (6.1) 
a 3 15 35 
Ym is the required mean radius. 


The external gravitational potential, V, expressed in spherical polar coordi- 
nates is 


fm f 
RR 


If this is re-written in terms of rm, we have 


M 3C- 
V = Lal ie =) © Pe(cos 6)+ .. , 


tm\R .R/ Mra? 


The ratio 3(C — A)/2Mrm2 will be called J*. 


5-(3) 


™m 
Using equation (4.1), 
I 2 I 8 172 
ta e(r + 20-4) —=m(1 + e+ Het), (6.2) 
\ 6 15 2 21 735 


To obtain the gravity formula, equation (3.7) is, with the use of equation (2.2), 
written as 


8 
(: ‘cae 2) ++m(5+—e+ Ba) costs | x 
ge eer SN 7 49 


x E - ~ (2+ 3e + 4e) cos?d + 5 e%(1 + 3e) cos#} — 28 cost |. 


By expressing this in terms of Legendre functions of cos $ and taking the con- 
stant parts, the value of gm is obtained: 


{M 2 I I 7 22 
gn = (1 -=m—*e-*e- Zam + =), 
ab 3 3 5 20 105 


or, in terms of the mean radius, 


see én). (6.3) 
135 $40 


It is an advantage that this relation does not involve e but only higher powers. 
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The complete gravity formula is 
20 I Pa. wb 
27 94s 
4 982 


12 
x Po(cos 9) + |= — -—em + ——e*m —- — 
35 7 1155 


8 
+ ——(3e2m — 5e%)P¢(cos 9). 
231 


7. Summary 
The terms of order e° in the gravity formula derived in section 3 are 
sel m( + aid m2 — Sn ~ . =") sin’d + — em(— e “ m)sin?ag 
8 4 98 8 7 2 
I 


15 . : 
~ e2/ —m— e| sin?d sin?24 ‘ 
8 \2 
With e = 0-003 37 and m = 0-003 45, the coefficients of the terms are respec- 
tively: 
of sin?¢d +0°24 mgal; 
of sin?2¢ — 0-02 mgal; 
of sin?¢ sin?2¢ + 0-03 mgal. 
It is improbable that the external field will be analysed to this accuracy for 
some time to come. 
The formula for J is 


I I 2 11 
e(: ~ -e} - —m( 1 —-e+ —¢). 
se... 7 49 
The relative magnitude of the final term is 


. oo 
2 49 
that is about 2-5 x 10-6. 
The relative precision that could be achieved at present in the determination 
of J from satellite observations is about 1 x 1074. 


13 Percy Road, 
Hampton, Middlesex. 
1959 March 16. 
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Tidal Friction from Currents near the Seabed 


H. Charnock 
(Received 1959 May 6) 


Summary 


The current in the lowest two metres of a tidal stream has been 
found to increase linearly with the logarithm of the height above the 
bottom. Independent measurements of the shearing stress on the sea- 
bed allow a determination of von Karman’s constant which was found 
to be approximately 0-4. 


1. Introduction 


Estimates of the shearing stress of a tidal current on the seabed have been 
published by Bowden & Fairbairn (1952, 1956). The earlier paper represents 
the outcome of attempts made over many years to determine the drag on the 
bottom using a dynamical method. It is made clear that this method is technically 


difficult and time-consuming and that the results are difficult to interpret. The 
observations described in this paper were made to supplement their investiga- 
tion and to establish a simpler but less direct method. 

From observations of shear flow in the laboratory, it has been confirmed that 


the velocity profile close, but not too close, to a boundary can be represented by 
the well-known expression 


U z 2 

a7 yee (1) 
where U is the velocity at a height z from the boundary; U, is defined by 
to = pU,2 where 79 is the shearing stress on the boundary and p is the density of 
the fluid; k, von Karman’s constant, is about 0-4. For flow over a rough boundary 
zo, the roughness parameter, depends only on the geometry of the surface. Over 
a smooth boundary the effective zo is proportional to v/U, where v is the kinematic 
viscosity. 

This expression, with the same value for von Karman’s constant, has also 
been shown to represent the vertical profile of the mean wind near the ground in 
conditions of neutral hydrostatic stability, so that it seemed likely that it would 
also represent the flow near the bottom in a tidal stream. Logarithmic profiles 
have been described by Lesser (1951); the present observations were made in 
the area where Bowden & Fairbairn had already made independent observations 
of the shearing stress, in the hope of confirming that the laboratory value of von 
Karman’s constant was again applicable. 

There are advantages of working in shallow-water tidal currents; they provide 
an almost neutrally stable flow of high Reynolds number in which the shearing 
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stress varies almost linearly from the bottom to the top. There is also a relatively 
long period variation of the mean velocity with time. Against this one must set 
the disadvantages of working from a ship; robust instruments are necessary and 
when in use neither they nor the lower boundary can be inspected. 

The results have practical applications not only to tidal dynamics but to the 


increasing number of problems involving the off-shore disposal of radioactive and 
other waste. 


2. Apparatus 


An instrument for measuring the current at various distances from the bottom 
was constructed and its general arrangement can be seen in Figure 1. It consists 


Fic. 1.—Outline drawing of current measuring stand. The vertical 
tube is 3°8cm in diameter. 


of up to five cupwheels centrally mounted at heights up to 2m on a vertical steel 
tube which stands on a tripod base. The cupwheels each have three aluminium 
alloy cups of almost conical shape and with beaded edges. They rotate on simple 
ball bearings in Tufnol housings. Stainless steel balls and Monel races resist the 
seawater corrosion fairly well. 

Each cupwheel carries two small magnets by which a simple waterproof switch 
(see Figure 2) is closed twice per revolution. The magnetic switches are con- 
nected by electric cable to a bank of telephone message registers on board the 
anchored ship. 

Some difficulty was experienced initially in preventing the electric cable and 
the lowering wire from fouling the cups when the ship rolled. The method finally 
adopted was to lower the stand on 4mm diameter galvanised steel wire rope to 
which was attached a small pellet or trawl float a few feet above the top of the 
stand. The float supports the lowering wire above the cupwheels. The electrical 
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cable was led from a junction box along one of the legs of the tripod and brought 
to the ship in a large bight. 

The instrument was calibrated on a whirling arm at the Admiralty Research 


Laboratory, Teddington and in the water channel in the Fluid Dynamics Depart- 
ment of Liverpool University. 


3. Observations 


Preliminary observations were made from the m.v. William Herdman in 1955 
July. The site was Red Wharf Bay, an area off Anglesea where Bowden & Fair- 
bairn (1952) had made estimates of the drag coefficient of the seabed. A full 
description of the area can be found in their paper. Supporting observations 
were made in the Menai Strait which was used as a shelter from bad weather. 

















Fic. 2.—Outline drawing of magnetic switch. 
A, two-core cable sealed into brass plug B; C, O-ring seal; D, soft iron 
block; E, silver contacts; F, tubular brass casing, 1-27 cm in diameter; 
G, insulating pad. 


The uppermost cupwheel of the four in use was lost at an early stage so that 
most of the profiles were based on only three points of which the lowest did not 
give consistent results. ‘The observations are shown in Figures 3a and 3b in which 
the mean (10 min average) velocities at 199, 62 and 31-5 cm are plotted against the 
velocity at 1oocm. There is an indication that the velocities are linearly related 
and their ratios are plotted against the logarithm of the height in Figure 3c. Apart 
from the bottom point, which has been ignored in drawing the line, the profile is 
roughly logarithmic and the roughness parameter appears to be about o-1cm. If 
k in Equation 1 is taken as 0-4 then the corresponding value of the drag coefficient 
based on the 100cm current (Cjo9 = (U,./Ui00)), is 3-4 x 107%. 

In view of the instrumental deficiencies it was necessary to repeat the observa- 
tions and an opportunity arose in 1957 July. The measurements were made in the 
Red Wharf Bay area, again from the m.v. William Herdman. On this occasion, 
thanks to improved technique, all five cupwheels in use survived and individual 
profiles could be drawn. These are shown in Figure 4 in which the mean velocity 
is plotted against the logarithm of the height above the bottom. Most of the profiles 
are approximately logarithmic; the occasional discrepancies could possibly be 
explained as due to turbulent fluctuations but are more likely to be due to jellyfish 
which were regrettably numerous on this occasion. The lines drawn on Figure 4 
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Fic. 3(a) Plot of Usi-s, Use, Uiss against Uioo, Menai Strait 1955 July 2s. 
(b) Plot of Usi-5, Use against Uioo, Red Wharf Bay, 1955 July 26. 
(c) Plot of [Uz/Ui00] against log z. 

x—Menai Strait observations. @—Red Wharf Bay observations. 
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Fic. 4.—Current profiles, Red Wharf Bay, 1957 July s. 
The profiles a, b . . . | were observed at successive times. They are 
plotted on the same scale but relative to arbitrary origins, the values 
of U2co being given in Table 1. 
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have been drawn by eye and values of U,, zo, etc. deduced from them using 
k = 0-4 are given in Table 1. Values computed using the least square technique 
do not differ significantly. It appears that the roughness parameter was about 
o-3cm and this is confirmed by the profiles of mean velocity ratio illustrated in 


Figures 5a and b. The corresponding value of the drag coefficient based on the 
100cm current is Cjo9 = 4°7 x 10°73. 
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Fic. 5(a) Plot of Uso, Usei, Uzs, Uisz7 against U200, Red Wharf Bay, 
1957 July s. 
(b) Plot of [Uz/U200] against log z, Red Wharf Bay, 1957 July 5s. 


Discussion 


The large amount of long period large-scale turbulence shown by the time 
variation of Usoo (‘Table 1) is noteworthy. 

The flow is hydrodynamically rough and the value of 0-1 to 0-3cm for the 
roughness parameter (Cjo9 = 3°4x 107% to 4°7x 107%) is longer than would be 
expected for smooth sand; it is thought that the bottom was roughened, perhaps 
by sand ripples. 

To compare these results with those of Bowden & Fairbairn (1952) it is neces- 
sary to transform their value for the drag coefficient to make it comparable with 
Cioo. Their analysis relates the amplitude of the semi-diurnal component of the 
tidal velocity to the amplitude of the semi-diurnal component of the frictional 
force; the implied linearization can be regarded as replacing |sin ot| sin ot by 
(8/37) sin ot and neglecting higher harmonics (see Proudman 1953, p. 309). 
Their value for the drag coefficient (1-8 x 10) should therefore be multiplied by 
37/8 to make it comparable with our measurements. 

A further adjustment is necessary because their drag coefficient is based on 
the mean velocity over the whole depth and not on the velocity at roocm. They 
give the mean ratio of the 100cm and depth mean velocity as 0-77 so that their 
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estimated value for Cioo is given by 
1°8 x 10°3 x 37 
(0°77)? 8 


If all these observations are taken at their face value the value of k was 0°39 
on the earlier and 0-46 on the later occasion. It must be remembered that all 
three sets of observations are subject to observational error, were made at dif- 
ferent times and not in precisely the same place, although the bottom is fairly 
uniform in the area. It seems reasonable to conclude that the laboratory value of 
von Karman’s constant is applicable in tidal streams and that observations of 
near-bottom currents can provide acceptable estimates of tidal friction. 

In a later paper Bowden & Fairbairn (1956) made estimates of the shearing 
stress on the bottom of Red Wharf Bay from the covariance of vertical and hori- 
zontal fluctuations of current. These lead to estimates of drag coefficient which 
are about half those obtained above. It is possible that some of the momentum 
flux was missed by the instruments being of insufficiently rapid response or by 
the observing periods being too short. But a more likely reason for this underesti- 
mate is that the corresponding mean currents were not measured directly but were 
estimated from earlier work. ‘The instrument described in this paper is therefore 
a useful adjunct to the electrode flowmeter which they have used for their fluctua- 
tion measurements; simultaneous observations of mean flow and fluctuations 
have been made and are now being analysed. 





Cus = = 36x 1073. 
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Reports on the Progress of Geophysics 


Developments in Dynamical Geodesy* 


A. H. Cook 


Summary 


Progress in dynamical geodesy over the last ten years is reviewed 
with particular attention to the significance of observations of the 
orbits of artificial satellites and of direct measurement of the distance 
of the Moon by radar. As a result of these observations together with 
a great increase in the number of gravity observations, especially at 
sea, and with more extensive areas of geodetic triangulation, the 
ellipticity of the Earth can be determined to one part in 3000 and 
the radius to a few parts in 100000. It is suggested that the Earth’s 
gravitational field should be specified by the parameter J in the spherical 
harmonic expansion rather than by the derived ellipticity and that 
the mean rather than the equatorial values of the radius and sea-level 
gravity should be employed. 


1. Introduction 


Dynamical geodesy is the study of the overall shape and size of the Earth 
and particularly, of the properties of its gravitational field. Its history begins 
with Newton’s study of the motion of the Moon and the shape of the Earth and 
the problems have engaged many of the most able mathematicians. Until the 
beginning of this century, theory was far ahead of observation because of the 
problems of making measurements that were sufficiently exact (relative errors 
should be less than one part in a million) and were sufficiently numerous and 
evenly spaced over the Earth. The last few years have seen remarkable changes 
in this situation, of which some account is given in this review. 

The period covered is the ten years or so since Sir Harold Jeffreys published 
his general adjustment of all the data then existing concerning the shape and size 
of the Earth (Jeffreys 1948) and to appreciate the significance of the new develop- 
ments in this period it will be helpful to summarize that paper. The observations 
which Jeffreys used were: 

(a) Measurements of arc lengths and angular amplitudes of triangulation net- 

works—these give the size and ellipticity of the sea-level surface. 

(b) Free-air anomalies of gravity. 

(c) The luni-solar precession of the Earth. 

(d) Certain aspects of the motion of the Moon. 

(e) The absolute value of gravity. 

(f) The trigonometrical parallax of the Moon. 


* Communication from the National Physical Laboratory 
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Except for the Earth’s precession, these all involve observations relative to the 
external gravity field of the Earth and in particular, relative to the equipotential 
surface which, over the seas, coincides with mean sea-level. This surface will be 
referred to rather loosely as the geoid. The geometrical form of this surface is 
related to the external gravity field and, in particular, to the values of gravity on 
the surface. We are ignoring in this review the details of the local irregularities 
of the geoid and the gravity field so that these irregularities have the effect of errors 
in the determination of the gross parameters of the geoid and gravity field. 

In the measurements of arcs of triangulations the measurements of the arc 
lengths can easily be reduced to the geoid since the heights of the observed lines 
above the geoid can be found from spirit levelling or otherwise. The effect of 
irregularities in the geoid will then be to make the measured distance slightly 
greater than that for the best fitting spheroid. This is only a second-order effect. 
The astronomical observations required for the angular amplitude of the arc are 
made relative to normals to the geoid and irregularities in the geoid affect the 
results directly. These irregularities are correlated and Jeffreys discussed in some 
detail the allowance to be made for them in estimating the uncertainties of the 
results. 

The variation of gravity over an ellipsoidal equipotential surface is related to 
the geometry of the surface by a well-known theory (Darwin 1899, Jeffreys 1952, 
Pizetti 1894) which may be extended to cover small irregularities of the surface. 
The geoid is very closely an exact spheroid but the topographical surface on which 
gravity observations are made is not an equipotential surface. However, it is known 
that to a good approximation, the form of the geoid may be found from free-air 
anomalies. In earlier papers, Jeffreys (1941b, 1943) had studied the spherical 
harmonic analysis of free-air anomalies and he incorporated the results in his general 
adjustment. ‘The parameters he used in that adjustment thus depend essentially on 
observations of gravity prior to 1939 and the distribution of these observations was 
not very satisfactory for estimating the amplitudes of the various spherical harmonic 
components since large gaps in the data meant that the estimates of the various 
coefficients were not independent. 

The Moon moves in the Earth’s gravitational field and certain periodic perturba- 
tions in her motion are due to the second harmonic variation of the Earth’s gravity. 
The quantities are however, very small and the effects depend also on the Sun’s 
gravitational attraction and on the figure of the Moon. In consequence, the stan- 
dard deviation of the second harmonic in the Earth’s field as found from lunar 
observations is about 5 times that found from surface gravity observations. 

The luni-solar precession of the Earth is determined essentially by the ratio 
(C — A)/C (Cis the moment of inertia of the Earth about its polar axis, A that about 
an axis in the equatorial plane) and the mass of the Moon, although there are also 
small terms arising from the planets and from general relativity. The coefficient 
of the second harmonic in the external gravitational potential is proportional to 
(C—A)/Ma?, where M is the mass of the Earth and a the equatorial radius, what- 
ever the internal constitution of the Earth, but if the stresses within the Earth are 
exactly hydrostatic, then the ellipticity of the surface may be inferred from the 
ratio (C—A)/C (Darwin 1899, de Sitter 1915b, Bullard 1948, Jeffreys 1952). 
At the time when Jeffreys made his analysis, the uncertainty in the value of 
the ellipticity so derived, due to departures from the hydrostatic state, was com- 
parable with the observational uncertainties of other data from which the ellipticity 
may be derived. 
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The absolute value of gravity does not enter the determination of the ellip- 
ticity from the variation of gravity with latitude because all the information that 
can be obtained from that variation is contained in the ratio of gravity at one 
latitude to that at another and relative measurements of gravity, particularly with 
pendulums, yield just such ratios. The absolute value of gravity is however re- 
quired in determining the size of the Earth from the lunar parallax. The value of 
the Earth’s gravitational attraction at the Moon’s mean distance R is equal to 
Rw? where w is the Moon’s mean angular velocity. The ratio of this quantity to 
the mean value of gravity gm at the surface of the Earth is, by the inverse square 
law, %m?/R2 where rm is the mean radius of the Earth. 

But the Moon’s mean horizontal parallax, 7’¢, is 2r7m/R and so 


2w2 
‘Ss 


(These expressions are only approximate, see Section 7.2). 

The main uncertainty in applying this equation arises from the measurement of 
the Moon’s parallax. This depends on observations of the Moon at Greenwich 
and the Cape Observatories and is affected by local deflexions of the vertical at 
these places as well as by the uncertainty in the distance between the observatories 
in terms of the Earth’s mean radius. The agreement between the lunar parallax 
as directly observed and as calculated dynamically from survey measurements of 
the Earth’s size has in consequence not been too satisfactory. 

By paying careful attention to errors of observation, particularly those which are 
not independent, Jeffreys was able to derive values of the Earth’s equatorial radius 
and ellipticity which satisfi d almost all the data. 

His final values were: 


a = 6378-099 + 0'116km, 
el = 297°10+ 0°36, 


the uncertainties being standard deviations. 

The largest change from previous work was the reduction of a by 0:28gkm 
from Hayford’s (1910) value, due to the use of free-air instead of isostatic anomalies 
for the reduction of the astronomical observations. 

It is now clear that Jeffreys’s values cannot.be accepted any longer. This is 
the consequence of three major technical advances made in the years since 1945. 
The first of these is the development of methods of measuring gravity rapidly 
and accurately on land by spring-type gravity meters and at sea by improved 
submarine pendulum apparatus and more recently by ship-borne gravity meters. 
Although such methods had given results before 1940, they have been so intensively 
developed and exploited since 1945 that knowledge of the world-wide features of 
the Earth’s gravity field has increased many-fold. 

The second development is the launching of artificial satellites. The motions 
of these bodies provide estimates of the Earth’s ellipticity and other parameters 
of the external gravity field that are more accurate than Jeffreys’s estimates from 
surface gravity, even though so far only three satellites, and those not with the most 
favourable orbits, have been observed with the necessary precision. 

Finally, it seems that successful radar observations can now be made of the 


Moon, giving her distance with much greater accuracy than can be obtained from 
the visual parallax. 
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In addition the classical method of arc measurements has been extended to 
much greater areas than heretofore with observations of higher accuracy. 

This review is devoted principally to the effects of these major developments 

but recent work in other fields is also considered. It is first necessary to discuss the 

theoretical basis of the study of the Earth’s gravity field and this is done in the 


next section, attention being paid particularly to the meaning and definition of the 
parameters involved. 


2. Theory of the Earth’s external gravity field 
The theory of the Earth’s external field comprises the following topics: 


(a) Expansion of the potential of an arbitrary body in spherical harmonics. 
This gives the form of the potential required in the theory of satellite 
orbits. 

(b) Figure of a rotating body with specified moments of inertia bounded by a 
spheroidal equipotential surface. 

This problem admits of an exact treatment and enables the ellipticity of the 
surface and the variation of gravity over it to be related to the moments of inertia 
and the external potential. 

(c) Distortions of the geoid corresponding to higher harmonics in the potential. 


(d) Application of (b) and (c) to the actual Earth which is not bounded by an 
equipotential surface. 


2.1 Expansion of the potential in spherical harmonics 
The basis of this expansion is MacCullagh’s theorem, according to which the 
first terms of the external gravitational potential V of an arbitrary body are 
{M (A+B+C-3I) 
=—— + + oe 
R f 2R% 
where fis the Newtonian constant, 
M is the mass of the body, 
A, B, C are the principal moments of inertia, 
R is the distance from the centre of mass to the external point, 


J is the moment of inertia about the line joining the centre of mass to the 
external points. 


V 





For the Earth, which is highly symmetrical about the polar axis, C is taken 
to be the moment of inertia about that axis and A and B to be the (nearly equal) 
moments about the principal axes in the equatorial plane. 

Then if @ is the co-latitude of the external point, 


I = Asin20+C cos?6 
and 


{Mitm food iain 
$5 Tn R R Mrn?2 

provided R exceeds the greatest radius vector of the surface of the body. 
Tm is the mean radius of the Earth, the definition of which is discussed below. 
It is usual to express V in terms of (a/R) where a is the equatorial radius but there 
are good reasons for using 7m instead, as was first indicated by de Sitter (19152). 
Most observations relating to the size of the Earth—geodetic measurements of 


Po(cos 8) +... ), 
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arcs, visual lunar parallax, radar observations of the Moon—are made in middle 
latitudes and the quantity which is most directly involved is either the radius of 
curvature of the arc or the radius vector of the observatory. It is desirable to state 
the theory in terms of a parameter as nearly equal as possible to the distance con- 
cerned and not one, such as the equatorial radius, which is derived from them by a 
relation involving the ellipticity. Furthermore, the quantity m, equal to the ratio 
of centrifugal force at the equator to gravitational attraction, is almost independent 
of e when expressed in terms of mean gravity and mean radius. Formally it does 
not matter what parameter, @, 7m or some other, is used in MacCullagh’s theorem, 
although a change in the definition of J (that is }(C — A)/Ma?) is also involved, and 
in this review the theorems will be stated in terms of rm. 


2.2 Figure of a rotating spheroid 


The theory of the figure of a spheroidal equipotential surface bounding a 
rotating body has been established by two methods. In the first (Darwin 1899, 
de Graaff Hunter 1935) the external potential has been expanded in a series of 
spherical harmonics. Is it open to the criticism that spherical harmonics are not 
orthogonal in the space outside a spheroid but as far as terms of order e? the same 
results are obtained by the second method (Pizetti 1894) in which ellipsoidal func- 
tions orthogonal outside an ellipsoid are employed. Alternatively, spheroidal 
harmonics may be used, as has been done in carrying the theory to include terms 
of order e* (Cook 1959), an extension that seems desirable in view of the latent 
precision of satellite observations. 

The important results are the following; they are stated in terms of the mean 
radius and mean gravity defined as 


rm = a(1— }e— }e*) 
and 


_fM{ 2 2 / 19e ) 


2m po — oe f 
;. BEG 
These are the mean values of radius vector and gravity expressed as functions 
of geocentric co-latitude. 


Gravity formula in terms of ellipticity: 


9 
— 


Tm 


5 15m2 17em\ . I : 
&lge = 1+ (=m—e +——- ) sin? + —e(e— 5m) sin22¢; 
2 4 14 | 8 
I 8 I 2 
glem = ttne+ m+ —Sem— nt (2m—e— eae m+ >m?) sin’d + 
oe, Se ee 2 3 4 6 


I 
+ gue 5m)sin22¢. 
Relation between ellipticity and moments of inertia: 


We TRS 
J = e-—-m-—-e? +-em; 
ake 


j* = e——m +e em; 
2 21 
J®+—-m—— J*2 4m? +—J®m. 
Ss | 56 389 42 
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Gravity formula in terms of moments of inertia: 


3 6 63 6 


Em 


I II 283 ; I 9 
+(2m—]* ——]** ——]*m+ —-m%) sin?¢ + J *? — 4.]*m —-m?) sin®2¢. 
6 7 56 8 4 


The symbols are as follows: 


¢ is geographic latitude, 


(w is the angular velocity of the Earth) 
ge is the value of gravity at the equator, 
Be a RP 


= 1+-e+=—m+-e? —-—em—-=m?, 


8m = 2. ee. oo 
2.3 Application to the actual Earth 


Whereas a satellite moves in the potential field of the actual Earth with an 
irregular surface, so that from its motion the true moments of inertia of the Earth 
are obtained, the formula with which surface gravity is compared applies to an 
idealized body bounded by an equipotential surface. To first order, the actual 
gravity values may be brought into correspondence with the idealized ones by 
applying the free-air correction—it has been pointed out that the significant 
quantity here is the potential difference between the point of observation and the 
geoid and not the gravity difference (Cook 1953). The more exact theory of this 
correspondence probably requires to be re-examined especially as a second 
harmonic in the distribution of land might introduce an error into the estimated 
second harmonic of free-air anomalies. The usual approach is to set up a model 
Earth with the correct second-harmonic field but with the higher harmonics giving 
the correct value of gravity at the external surface only. This avoids a number of 
difficult questions and is sufficient for a first-order treatment of the departures 
from the second harmonic field. 


3. Spherical harmonic analysis of gravity values 


The principal difficulty in the analysis of gravity values is the uneven distri- 
bution of observations together with the correlation between gravity values at 
places separated by very considerable distances, 

The net effect is that the coefficient of any one harmonic depends to some 
extent on all the others whereas if the observations were uniformly distributed all 
the coefficients would be independent. The situation has greatly improved since 
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Jeffreys made his analysis (1941b). He then had available fewer than 10 000 obser- 
vations and almost the whole of the South Pacific and the southern oceans as well 
as the polar regions were unsurveyed. Nowadays many land areas are covered 
quite densely with observations made with gravity meters, some seas especially the 
Mediterranean, are adequately surveyed and the observations in the North Atlantic 
and the North Pacific oceans have been greatly extended. Observations have also 
been made in Canada as far north as latitude 82°30’. The main gaps now lie in 
the southern oceans and Australia; in addition, Russian results are not generally 
available*. Very many observations indeed have been made in some land areas 
and the problems of handling all the data and reducing them to forms suitable 
for geodetic use are very considerable. One would like to have mean free-air 
anomalies in areas one degree of latitude by one degree of longitude, together with 
the mean heights of these areas, for it is necessary to analyse the variation of free- 
air anomaly with height. Two groups are working on this problem—the Inter- 
national Gravity Bureau in Paris which, until his untimely death was under the 
direction of Dr Lejay, and Professor W. A. Heiskanen’s group at Columbus, 
Ohio. Few results have so far been published although Professor Heiskanen has 
summarized his work (Heiskanen 1957). Even now, only about 15 per cent of all 
one-degree squares have observations. 

The effect of correlation and the uneven distribution of observations was 
thoroughly discussed by Jeffreys (1941b, 1943) while Kaula (1957) has developed 
an alternative approach which may be applied to the problem of spherical harmonic 
analysis. The results of analyses are expressed in the form 


g = ge1+f sin*d— 5-9 x 10-6 sin22¢4). 
The term with sin22¢ is calculated theoretically from the condition that the 
equipotential surfaces are exact spheroids. g, and 8 are determined from the 
observations. Some results are given in Table 1. 
Table 1 
Results of spherical harmonic analysis of free-air anomalies 


Reference &e B e} e 


Jeffreys, 1948 978-o51cm/s* 52863 x 107° 296-85 1-640 x 10° 
+0°0075 +0-66 


Heiskanen, 1957 978-0497 5°2902 297°2 1°633 
Krassovski, 1940 978-049 5°3029 298-4 1°622 


Zhongolovitch, 1952 978:0573 §°2837 296°61 1642 
+0:0064 0°57 


Heiskanen’s result has been criticized by Jeffreys (1959) for neglect of the 
effects of correlation. 

There are two theoretical points in the analysis of gravity observations that 
require study. One is the effect of low order harmonics in the height of the topo- 
graphy which may come from two causes, namely from any error in the free-air 
correction formula and also, through the imperfect formation of means, from the 


* But see Zhongolovitch 1952. 
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correlation of free-air anomaly with height. The second point is that anomalies 
are always given in terms of geographic latitude. However, since they are referred 
to a surface which is a spheroid they should be analysed in terms of functions that 
are orthogonal on a spheroid—that is, Legendre functions of confocal coordinates 
(see Cook 1959). 


4. Motion of an artificial satellite 
4.1 Theory of perturbations 


Within recent years a great deal of work has been done on the motion of a 
satellite in the Earth’s gravity field. Some of the earlier studies are too restricted 
or are in error and present determinations of e and J from satellite orbits depend 
on the work of King-Hele (1958), of Herget (1957) and of Krause (1956). 

King-Hele has carried the theory as far as terms of the order of J?. He takes the 
gravitational potential of the Earth in the form 


{M 2 8 
—{a/r —(a/r)8 —] Po(cos 6) — (a/r)® —D P,(cos 6)}, 
a 3 35 


where r is the radius vector and @ the co-latitude of the satellite, and a is the 
Earth’s equatorial radius. He then finds by a perturbation method that the rate 
at which the plane of the orbit rotates about the polar axis of the Earth is given by: 


[J+ J%all)*(— sin?a— 1) — GD/14) (a/P)°(7 sin®a — 4) + OU*e)]; 


the rotation is in the opposite sense to the motion of the satellite in its orbit. 
Here R is the semi-major axis of the orbit, / the semi-latus rectum, e the eccentricity 
of the orbit and « the mean inclination of the orbital plane to the equator. 

The other principal secular effect is the rotation of the major axis which 
King-Hele finds to be given by 


d, 
& = Key'2)Mald%a\R)' (5 costa 1)+O(Je) +OU2) 
where y = {M/a?. 

When cos®« = 1/5, the rate of rotation is zero; for equatorial orbits, the axis 
rotates in the same direction as the satellite and for polar orbits in the opposite 
direction. 

King-Hele considered just the second and fourth harmonic terms in the Earth’s 
potential but results of the spherical harmonic analysis of free-air anomalies do 
not exclude the existence of other harmonics. 

R. H. Merson and R. N. A. Plimmer of the Royal Aircraft Establishment are 
extending the theory, by a method somewhat different from King-Hele’s, to take 
account of other harmonics and find, for instance, that if there is a third harmonic 
with coefficient J3, it gives rise to a component of the rotation equal to 


— 3(f M/R*)*(a/] Jse sin w cot «(15 sin®« — 4). 
(Merson & King-Hele 1958). w is the mean argument of perigee. This rotation is 
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evidently very small for nearly circular orbits so that J3 can only be well deter- 
mined from highly elliptical orbits. J. de Graaff-Hunter is also studying orbital 
theory. American workers have used the theory developed by Krause (1956) 
and that of Herget (1957). Krause neglects the J? term yet keeps the D term which 
is of the same order. 

Superposed on the effects of the Earth’s gravitational field are those due to 
the resistance of the air. ‘The main effects of the atmosphere are to reduce the major 
axis, to reduce the eccentricity of the orbit and to reduce the period. These are 
quite different from the effects of the oblateness of the Earth and are readily 
separated from it by observation. The theory of the effects of air drag has been 
studied by many authors and so far there is no indication that drag will cause 
rotation of the plane of the orbit or of the major axis. A lateral force on the satellite 
due to rotation of the atmosphere might produce a rotation of the orbital plane and 
this would be revealed most directly by an orbit passing over the poles; another 
way of testing for such an effect would be to have satellites going westwards as 
well as eastwards round the Earth. (King-Hele & Merson 1958). 

So far the attractions of the Sun and Moon have been ignored; they cause 
perturbations amounting to about 1 in 10000 of those due to the oblateness of the 
Earth and it may eventually be necessary to allow for them, particularly if orbits 
with appreciably greater axes are used. There is also a very small effect from 


general relativity. 
4.2 Reduction of satellite observations 


A term in the gravitational potential that is periodic in longitude will not pro- 
duce any secular effect on the satellite orbit unless the period happens to 
coincide with that of the satellite. Such effects are ignored at present and we con- 


sider only the perturbations due to terms in the potential proportional to the 
Legendre functions P2(cos @), P3(cos @)... with coefficients Jo, Js, ...*. Then 
in general, the orbital plane will rotate at a rate 


C2 Jo+C3Ja+Cs Jar... 


while the major axis will rotate at a rate 


do Jo+d4 Jat... 


where the coefficients C2... dz... vary with the parameters of the orbit and in 
particular, C4 vanishes for sin’ = 4/7 and d2 for cos’« = 1/5. 

In favourable circumstances observations of both quantities for one satellite 
might enable J2 and J4 to be found separately. More generally, it should be pos- 
sible to determine a number of the J’s by combining the observations on different 
satellites. At present it is easier to get accurate observations of the rate of rotation 
of the orbital plane than that of the major axis, since the plane is the better defined. 

The Russian satellite, Sputnik 2 (19578), has been observed optically by kine- 
theodolite (Merson & King-Hele 1958) and the rate of rotation of the orbit is 
quite well defined. The method of analysis is described by King-Hele & Merson 
(1958). The result can be expressed as follows: 


1-0008j —0°312d = (— 6-8 + 0-3) x 10° 


where j = 5J/J, d = 8D/J and 8] and 8D are the differences from Jeffreys’s 
values J = 1-637 x 10-3 and D = 10-6x 107, 


* Note that J = —4 Je, and D = 3% J, in the conventional expression for the potential. 
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Sputnik 2 was also observed visually and photographically at eleven observa- 
tories in Czechoslovakia, mostly by plotting the position on a star chart. Buchar 
(1958) has determined the motion of the ascending node from these data and 
has calculated J from the result using a formula which agrees with King-Hele’s 
as to the terms proportional to J and D but from which the /? term has been 
omitted. In his numerical work, Buchar assumes that D has the theoretical value 
for a perfect spheroid, namely 10-7 x 10-6, and then gets 


J = (16278 + 0°0020) x 10-3 


corresponding to e~! = 297-9 + 0°18. 

The rate of rotation of perigee (major axis) of the orbit of Sputnik 2 is very 
small because 5 cos*« is close to 1. ‘The American satellites have a larger rate of 
rotation of perigee and J and D can be found separately by combining results from 
observations of this motion with observations of the motion of the node (rotation 
of the orbital plane). American satellites are observed by the Minitrack stations. 


Harvard announcement card 1408 of 1958 June 14 gives the following values of 
el, 


Satellite Observation ei 


19580 Node 298-0+0°3 
Perigee 297°8 

19588 Node 298-38 + 0:07 
Perigee 298-3 


These values appear to be derived on the assumptions that J? is negligible 
and the D has the value 10-6 x 10-6 for a truly spheroidal Earth. 
Cook (1958) has combined King-Hele and Merson’s result with Jeffreys’s 


(1943) values of the coefficients of the second and fourth harmonics in gravity, 
obtaining 


J = (1°628 + 0004) x 10-3 
D = (15-7 + 8-2) x 10-6 
This is not really a satisfactory procedure since it combines two sets of data 
that one would like to compare separately and better results are obtained by 


combining the results for Sputnik 2 and the American satellites. In this way 
Merson finds the following equations: 


From Sputnik 2 
J-—0o311D = 1°6227 x 1073 
from American satellites 
J +0°2645D = 1-6264 x 107% 
J+0°2035D = 16278 x 1073; 
the least squares solution is 
J = (1°625 + 0°001) x 1073 
D = (8+3)x 10% 
el = 298-2+0'1. 
* There is a misprint in the original paper. 





(Geophysical Discussion of the Royal Astronomical Society, 1959 January 23—see 
Geophys. J. 2 (1959) p. 251). 

Better equations for Vanguard 1 (Jacchia 1958) and Explorer 4 are now avail- 
able and the definitive results are those of King-Hele & Merson (1959) who have 
combined observations of the rates of rotation of the orbital planes of Sputnik 2 
(19578), Vanguard 1 (195882) and Explorer 4 (1958¢) and find: 


108J = 1624:6+0°3, 10°D = 6+1, 
10°8Jg = o-1r+1°5 and e-! = 298-20+0-03. 


To summarize the results so far obtained, the value of J seems to be deter- 
mined at about 1-625 x 107-3, D appears to be less than its conventional spheroidal 
value of 10-6 x 10-6 and J¢ is scarcely significant. 

In addition, observations on. the American satellite Vanguard 1 suggest that 
the term proportional to P3(cos @) in surface gravity has an amplitude of about 
5 mgal (O’Keefe, Eckels & Squires 1959). 

It is desirable that the results of satellite observations should be given primarily 
as coefficients of spherical harmonics in the potential rather than as the correspond- 
ing ellipticity because although the latter is clearly defined as (a—6)/a for a true 
spheroid, the definition can be ambiguous if the equipotential surfaces depart 
from spheroids. The coefficients of the spherical harmonics on the other hand, 


are unambiguous physical properties of the Earth and are the quantities that directly 
determine the motion of a satellite. 


4.3 Comparison with surface gravity observations 


Jeffreys’s analysis of free-air gravity anomalies gives 
103] = 1-640+0°007; e-1 = 296-85 + 0°66 


while the satellite observations give 
lol] = 16246010003; 1 = 298-2 +003. 


Jeffreys (1959) has discussed possible reasons for this discrepancy but in view 
of the fact that smaller values for e have been obtained from the same data (see 
Section 3) it is not unlikely that, when the more recent gravity observations have 
been incorporated into the analysis, better agreement will be obtained. A discussion 
of the available material on the lines of Jeffreys’s analysis is very desirable. 


5. The moments of inertia of the Earth 


The effect of the oblateness of the Earth is that the attraction of the Sun 
exerts a couple tending to rotate the Earth about an axis in the plane of its orbit 


and perpendicular to the line joining the Sun and the Earth and in consequence. 
the Earth precesses at a mean rate equal to 


3 n2 C-—A 
2 c 


cos 6, 


where n’ is the Earth’s angular velocity in its orbit, 
w is the Earth’s angular velocity about its polar axis 
and @is the angle between the Earth’s polar axis and the pole of the ecliptic. 
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The Moon also contributes to the precessional couple and gives a precession 
equal to 
3C-An’ 1 P 
— —— — — cos, 
a C 
where n” is the Moon’s angular velocity in her orbit, 
2 is the ratio of the mass of the Earth to the mass of the Moon, 


6’ is the angle between the Earth’s polar axis and the pole of the Moon’s 
orbit. 


The total luni-solar precession is about 50” per year of which about 15” is 
due to the Sun. 

There are some minor complications in the theory. In the first place, allowance 
has to be made for the effect of the other planets. There is also a small correction 
entailed by general relativity. The net precession due to the Sun and Moon 
is §0°-39 per annum. 

The value of u is required in order to calculate (C — A)/C denoted by H and this, 


and the constant of precession P have been redetermined by Rabe (1950a, b) 
from observations of Eros as follows: 


P = 5493°°793 
po) = 81-375; 
these lead to 
H = 3:27559x 107%. 


Jeffreys’s (1948) value was 3:27237 x 1073, 

If it is assumed that the stresses in the Earth are hydrostatic, then J may be 
calculated from H (Jeffreys 1952). Since H is exceedingly well known, the un- 
certainty in the value of J found in this way comes entirely from ignorance of the 
true state of stress in the Earth. Jeffreys (1948) estimated the uncertainty of J 
to be 3 x 10-6 which makes this estimate of J less uncertain than those from any 
data other than satellite observation. The satellite value is however more precise 
so that it may be used in conjunction with H to compute the moments of inertia. 

We have 

H = 0°0032756; 
J = 0°001624%, 
C-A 


yo tte 00010831. 
Ma 


Cc 


Ma §(J/H) = 03306, 


A 
We 0°3295. 


If the Earth were in a hydrostatic state, e~! would be 297-10, the corresponding 
value of J being 1-637 x 10-8. This value of e~! is greater than Bullard’s (19 48) 
value because of the reduction in H. 

The comparison of the satellite value of J with that computed from the 
luni-solar precession demonstrates that the stresses in the Earth are not hydrostatic, 
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as is of course evident from the existence of mountains and ocean basins as well 
as from the large-scale departures from isostacy. 


6. The motion of the Moon 


Just as with an artificial satellite, the plane and the major axis of the Moon’s 
orbit rotate on account of the oblateness of the Earth. The theory was given in 
its most complete form by E. W. Brown (1908). At the distance of the Moon 
the higher harmonics have no effect but the rates of rotation depend also on the 
oblateness of the Moon and further, are much smaller than the rates of rotation 
due to the Sun and planets. Thus it is difficult to obtain accurate estimates of the 
parts that depend on the figures of the Earth and Moon while the solution for the 
Earth’s oblateness depends on that of the Moon. Jeffreys (1937, 1941a, 1948) 
discussed these problems and used lunar data in his general adjustment but now 
that the Earth’s oblateness is so well determined by satellite observations, it is 
possible to substitute the value of J in the lunar equations and re-examine the 
values of the lunar parameters. In addition to the quantity J the following para- 
meters referring to the Moon, are required: 

, 32C’-A’-B 3 B’-A’ 
J 2 2M’? ’ 2 M’a’?’ 


re 


C’— A’ ee, 


RET OS ee 


L’ = J'+4K’. 

C’ is the moment of inertia about the axis of rotation, A’ that about the axis 
pointing to the Earth and B’ about that tangential to the orbit. a’ is the equatorial 
radius. 

The following equations are then found for the motions of the node (rotation 
of plane) and of perigee (rotation of major axis): 


Perigee 3897] + 380L'’—1192K’ = 6"-411 + 0°039 per year; 
Node — 3648] —460L’ = —6"-179 + 0-039 per year. 


(the motions due to solar perturbations are 146427” for the perigee and — 69672” 
for the node). 


Taking J to be 1-625 x 107%, 


L’ = 5:46x 1074, 


= 1°28 x 1074, 
and so 


J’ = 42x 1074. 


The value of 8 is related to the inclination of the Moon’s axis and de Sitter 
(1915c) gives 
B(1 +.0°0047 f) = 0:0006286 + 0-0000015. 
f is equal to «/8, that is to 1— K’/L’. 


f may be determined from observations of the Moon’s libration in longitude 
though until recently the results were very discrepant. Recent discussions of new 
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observations appear to give more concordant results and Jeffreys (1957) finds 


B = 0:0006266 + 0-000 0022 
and 


f = 0673 +0°0018, 
while Habibulin (1958) obtains 
f = 0°604+0°03 
which entails 
B = 0-000 6279. 


Now the above values of K’ and L’ lead to f = 0-77 which appears to be 
significantly greater than either Jeffreys’s or Habibulin’s value. 
There is a further difficulty. We have 


3 C’-A’ 

2 M'a”? 

<a 
nic 2 M'a’ 


which with the above numerical values gives 


’ 


and thus 


= L’/B, 


gq = 0°87. 


This however, implies that the interior of the Moon is considerably less dense 
than the outer parts (if the density were uniform, g’ would be 0-6). It seems there- 
fore, that there are inconsistencies in lunar observations which are not removed 
by an improved value of /. 


7. The size of the Earth 
7. 


1 The radius from arcs of geodetic triangulation 


The radius of the Earth has for many centuries been found from measurements 
of the distances (arc-lengths) between places of which the geographical latitudes 
are known from astronomical observations. Among such measurements of histori- 
cal importance are the French measurements of arcs at the equator and in 
Lapland which showed that the Earth was, as Newton maintained, an oblate 
spheroid, and the measurements of the arc through France and Spain from which 
the ratio of the Metre to the Toise was obtained. If the observations were made 
on the surface of a spheroid the errors of the calculated radius would be deter- 
mined solely by the errors of observation. In practice, the observations are made 
at places which may be some distance above sea-level and are referred to surfaces— 
the equipotential surfaces passing through the points of observation—which are 
not true spheroids. The height above sea-level can be measured by spirit-level 
observations and corrections can be applied to reduce the observations to sea- 
level but the irregularities of the equipotential surfaces cause more trouble. If 
gravity surveys extended over the whole surface of the Earth, the deflexions of the 
vertical could be calculated at all points along an arc and the corresponding geo- 
graphical latitudes on a spheroid could be calculated from the true observed 
latitudes. At present, however, this is not possible and the most serious diffi- 
culties arise not from those parts of the deflexions due to local mass irregularities 
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but from those which maintain the same sign over a large part of an arc of obser- 
vation. The reduction of the observations taking these correlated errors into account 
was discussed by Jeffreys (1948). A further difficulty is that the observed arcs do 
not extend from pole to pole and so instead of giving a and e separately, they are 
best considered as giving a relation between a and e. Then if e is found from 
other data, a may be calculated. In this way, Jeffreys obtained: 


a = 6378-099 + 0°116km 
for 


e-! = 297°10. 


In this result some weight was given to an equation from the Moon’s parallax 
(see Section 7). Jeffreys’s result differs from the value due to Hayford (1910) 
adopted for the International Spheroid largely because Hayford applied isostatic 
corrections to the observed latitudes, a process which may introduce systematic 
errors. Jeffreys’s discussion of the observations available up to 1948 is the most 
thorough; the results up to 1950 have been summarized by Strasser (1957). 

These earlier results have been superseded with the completion of two major 
arcs of geodetic triangulation, that in North and South America and that from 
Scandinavia to South Africa where observations on the arc of the 30th meridian 
have at last been finished. Using these new data, Chovitz and Fischer (1956) 


have made two solutions; one with isostatic reductions and the other without. 
They assume that e-! = 297+ 1 and obtain 


a = 6378-240km + 100 (free-air) 
and 
a = 6378-285 km + 100 (isostatic). 


The uncertainties arise almost entirely from the uncertainty of e~!. If e~! is taken 
to be 298-2 + 0-03 (Section 4.2) the values become 


a = 6378-126km + 50 (free-air), 


a = 6378-165 km + 18 (isostatic). 


The reduction in the uncertainties is the consequence of the higher precision of 
e~1, 


The difference of the free-air value for e~! = 297 from Jeffreys’s value is due 
to the much longer arcs that Chovitz & Fischer used. 

Mrs Fischer (1959) has subsequently made a new estimate of a from observa- 
tions of the height of the geoid above the reference spheroid. ‘The material comes 
from large areas in N. America, Europe and the USSR as well as from smaller 
regions in Africa and S. America. The data are effectively of the same kind as 
measurements along arcs of triangulation but cover areas rather than lines and the 
results correspond to the “ free-air” values above. For the conventional value 
e~1 = 297 (exact) 

a = 6378-253 km 


and with the satellite value e~! = 298-2 + 0-03, a will be somewhat less, about 
6378:200km. No firm estimate of the uncertainty has been given. 


7.2 The radius from the distance of the Moon 
The distance of the Moon is connected with the radius of the Earth as follows: 


Apart from perturbations due to the Sun and other causes, the Moon’s orbit 
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about the Earth is an ellipse of which the semi-major axis R is given by 


n”2R3 = {M 


where n” is the Moon’s mean angular velocity in her orbit and M is the mass of the 
Earth. Now the product fM is obtained from gravity observations on the surface 
of the Earth; in fact for a spheroid; 


2 2e / 19€ ' 
n-™ (2 sn} 
a2 2 


(The value of gravity should be increased by about 1 in 10° for the attraction of 
the atmosphere which does not enter surface measurements.) Thus the observed 
value of R when combined with the very accurately known value of n” and with the 
value of gm, gives an equation for rm which is almost independent of e. 

Until very recently, R had to be found trigonometrically from the Moon’s 
parallax w{. The definitive geometrical value is derived from observations of 
Christie & Gill (1911) at the Greenwich and Cape Observatories. However, the 
distance between these places had not then been measured directly and so had 
to be calculated from the geographical latitudes using assumed values of a and e. 
In addition there may be a systematic error in the parallax due to the deflexions of 
the vertical at the observatories. O’ Keefe (1950) in the United States, has effected 
a much more satisfactory trigonometrical determination from photoelectric ob- 
servations of occultations made at places the relative positions of which could be 
found directly from the geodetic triangulation of the United States without any 
assumption of values of a and e. He finds 


R = 384410km. 


More recently the distance of the Moon has been measured successfully by 
radar at the Naval Research Laboratory, Washington D.C. (Yaplee, Bruton, Craig, 
& Roman 1958) and at the Royal Radar Establishment, Malvern (Hey & Hughes 
1958). 

Although the results are only preliminary they show that it should be possible 
to obtain R to almost one part in a million, provided that the theory of the Moon’s 
motion is adequate. Brown (1905) considers that his calculated values of the lunar 
parallax are correct to 0”-o1, that is, to about 3 in 10%, so that the present theory 
may not be quite good enough. This depends to some extent on how the observa- 
tions are made. The effects of the Sun on the Moon’s orbit are to increase its 
size and to introduce periodic perturbations with arguments such as the difference 
of the longitudes of the Sun and the Moon. If the observations were to extend over 
a sufficiently long time, the periodic perturbations could be determined observa- 
tionally but if the observations cover only a short time interval, then the coeffi- 
cients of the periodic terms must be calculated theoretically accurately enough to 
apply corrections for them. In any case, the constant term should be evaluated 
theoretically to about one part in a million of the radius vector. 

Radar observations and occultation observations actually determine the distance 
Dons between the observatory on the Earth and a point on the surface of the Moon. 
Yaplee and others and Hey and Hughes both find that the effective reflecting 
surface of the Moon for 10 cm radio-waves is a fairly small area normal to the line 
joining the observatory and the centre of the Moon. If the length of the radius 
vector of this part of the surface of the Moon is assumed to be known, then the 

Q 





238 A. H. Cook 


equivalent distance D from the observatory to the centre of the Moon may be 
found. Now let r be the length of the radius vector of the observatory and @ 
the angle between it and the line joining the observatory to the centre of the Moon. 
Hence 


D2 = R2+r2—2Rr cos@ 


where R is the distance of the Moon at the time of observation. 

R can be calculated in terms of R and r in terms of r», and e. If R is replaced 
by its expression in terms of gm, 7m, e, and m, we have an equation for D involving 
these quantities (the non-spherical shape of the Moon should also be allowed for). 
D can be determined to about 1 in 300000 for a single value of @ and if the obser- 
vations are extended over a long period the resulting combined equation should 
have a weight rather higher than that of equations derived from surface measure- 
ments on the Earth. In this equation m has an uncertainty of 6 x 1075. gm depends 
on measurements of the absolute value of gravity as well as on relative measure- 
ments; measurements now in progress (see Cook 1957) should give the former to 
1 in 108 while the latter depends on e. Moon radar observations should therefore 
give an equation between rm and e with relative errors of the order of 1 in 1o® 
as compared with the equations from measurements over the surface of the Earth 
which can scarcely have greater precision than 1 in 10°. 

Moon radar observations for this purpose have to be made at a wavelength of 
about 10 cm. The available power decreases at shorter wavelengths while at 
much longer wavelengths the echoes fade in a complicated way and the refractive 
index of the ionosphere introduces large errors. The difficulties do not occur at 
1ocm. In the two experiments so far made the pulses have been about 5s long 
and a resolution of about } km has been achieved. 


Summary 


Since Jeffreys’s 1948 discussion, quite new ways of finding a and e have been 
developed and in addition, much new material for the classical methods of gravity 
survey and arc measurements has been acquired. The new gravity material has 
not been fully discussed so far and radar observations of the Moon have not given 
final results but very good values of J and e come from satellite observations and 
of a from arc measurements and geoidal heights. These are 


103 J = 1°6246 + 0°0003 
e-1 = 298-20+ 0°03 
a = 6378-200km. 


The uncertainties in 7m and gm are of the order of e2(1 in 108) and if with further 
improvements in technique, uncertainties of the order of 1 in 108 can be achieved, 
questions of the definition of parameters such as gm, Tm and e will need to be re- 
examined. Satellite observations, it must be emphasized, give physical parameters 
such as J whereas the reduction of all other observations involves in some way the 
use of a best fitting spheroid which has no physical existence. 

The launching of artificial satellites is stimulating new work on the overall 
physical properties of the Earth and, in turn, improved values of the astronomical 
constants of the Earth, Moon and Sun will be very valuable when it becomes 
possible to launch space probes into orbits designed to make specific observations. 
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Note added in proof.—W. M. Kaula (Statistical and Harmonic Analysis of 
Gravity, U.S. Army Service Tech. Rep. No. 24, March 1959) has made an 
analysis of gravity anomalies assuming that they form a stationary Markov popu- 
lation. He does not derive the lower harmonics from the gravity data but takes 
them to be determined from satellite observations. 


Note added in proof.—Y. Kozai has estimated that 10/3 = — 2-20 +0-08. (The 
Earth’s gravitational potential derived from the motion of Satellite 195882 
Smithsonian Astrophysical Observatory, Research in Space Science, Special Report 
No. 22, p. 1. March 1959.) 
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Geophysics in Pakistan 


H. I. S. Thirlaway 


Introduction 


The Geophysical Institute at Quetta is directed by Mr S. N. Naqvi, Director 
of the Meteorological Service. It may be said to have been established on 1951 
December 26, when the first accurately timed seismograms were recorded, but 
the specially built laboratories and instrument rooms (the foundations of which 
were laid in October 1952) were not completely occupied until September 1953. 
The Institute is 5 600ft above sea-level at the base of a 10000ft peak 5 miles 
west of Quetta. Laboratory and instrument room space of 12 000 square feet, and 
5 000 square feet of library, workshop, service and store rooms, were designed and 
sited by UNESCO advisers and built by the Pakistan Government on seismic 
resistant principles enforced by law in Quetta. 

Plans for the Institute were initiated in 1950 by Mr Mohammad Aslam, former 
Director Meteorological Service, and by Dr Malcolm Adiseshiah, Chief of the 
Technical Assistance Department, UNESCO (now Assistant Director General). 
The intention was to create a centre in Pakistan for the geophysical study of the 
Earth and its atmosphere. Atmospheric physics, geodesy, geomagnetism, and 
seismology were to be the basic fields of activity; geodesy and the absolute field 
geomagnetic survey being executed by the Survey of Pakistan. The basic agree- 
ment between Pakistan and UNESCO was signed in the autumn of 1950. It pro- 
vided four consultants* in the basic fields, pilot and demonstration equipment to 
the value of $35 000 and 6 fellowships. Each year since 1951 UNESCO has 
continued to support the Institute with consultants, equipment and fellowships. 
The Meteorological Service built and equipped the Geophysical Institute at a 
cost of about $120000 and in addition has set up seismological observatories at 
Karachi, Lahore and Warsak Dam (near Peshawar) in West Pakistan, and Chitta- 
gong and Sylhet in East Pakistan. It now plans to build a student hostel at Quetta 
and a similar Institute at Chittagong to serve East Pakistan. 

From 1955 the Survey of Pakistan, directed by Mr M. N. A. Hashmie, ex- 
panded its co-operation with the United Nations through another Specialized 
Agency, UNTAA, to include photogrammetry and the manufacture of survey 
instruments. Assistance in the geodetic work was then transferred to UNTAA for 
administrative convenience. New laboratories for the Geodetic Institute of the 
Survey of Pakistan are expected to be built in Quetta within the next two years. 
The photogrammetric work of the Survey of Pakistan at Rawalpindi is based on an 
air photo survey of all West Pakistan (save the high Karakoram region) on a 


* E. M. Fournier d’Albe, G. Norgaard, H. I. S. Thirlaway and K. A. Wienert. With the 
exception of Thirlaway, all served up to 1954 when responsibility for all observatory work at the 
Institute was handed over to Pakistan. Thirlaway will leave Pakistan in 1960. 
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contact scale of 1/40000 carried out by the Photographic Survey Corporation of 
Canada under a Colombo Plan project. 

The Headquarter offices and laboratories of the Geological Survey, directed 
by Dr E. R. Gee are also in Quetta, and a bi-lateral project with the International 
Co-operation Administration of the United States provides consultants from the 
U.S. Geological Survey, scholarships and equipment to assist in the recent large 
expansion of the Survey. The Survey has a Geophysical Branch with a staff of 
6 geophysicists. 

The three Departments are well located geographically for the study of Earth 
physics, since Quetta is situated in the Baluchistan mountain arc at the centre of 
the re-entrant structures. 

The Geology Department of the University of the Panjab, Lahore, has recently 
added geophysics to its courses and started systematic field work. ‘Though the 
Department’s senior staff of four are all from UNESCO,* it is hoped that local 
staff will be available by the end of 1960. 

Lahore is also the headquarters for the Groundwater Organization of the 
Public Works Department. This group are applying geophysics to the serious 
water logging and salinity problem in the Panjab. A geophysicist of the U.S. 
Geological Survey* is associated with the programme. 

A serious shortage of well-qualified geologists and geophysicists in University 
and Government service is due partly to the demands of some half-dozen oil 
companies which are actively prospecting the Indus Basin and East Pakistan. At 
least four promising geophysicists have left the country to seek better financial 
and professional opportunities abroad. Before 1951 however there were no geo- 
physical centres, laboratories, nor geophysicists in Pakistan at all. This article 
summarizes the results in Earth physics since that time and introduces the 
geophysicists responsible. 


Seismology 


The central observatory was installed at the Geophysical Institute in 1953 by 
Mr Moiduddin Ahmad with the assistance of Messrs R. A. Shamsee, A. H. Naqvi 
and A. Rehman. Mr Moiduddin subsequently installed the Karachi station, and 
Mr A. Q. Khan those of Lahore, Warsak Dam (Peshawar), Chittagong (new site) 
and Sylhet. The first station at Chittagong was started by Dr M. A. Choudhury 
in 1955. Mr A. Q. Khan, the present seismologist, was responsible for calibrating 
all the stations, and introducing amplitudes and magnitudes into the consolidated 
monthly bulletin of seismological readings. 

Quetta, Karachi, Lahore and Chittagong are equipped with Sprengnether 
1s period sets. Warsak Dam has a single horizontal 1s period Sprengnether and 
a vertical Willmore component, and Sylhet a Wenner accelerometer. Quetta and 
Chittagong are also equipped with Willmore sets, and Quetta with a Wenner 
accelerometer, a single Milne-Shaw, and two 20s period Sprengnethers. There 
is need for a long period vertical to match the horizontal components at Quetta. 
The maximum drum speed used is 1 mm/s. 

The Sprengnether sets have been found very robust, readily adjusted and 
easily maintained. The Willmore sets, though better for near earthquake work, 
have proved more difficult to maintain at peak efficiencv. They have been used 


* K. Helbig from Munich is the geophysicist. 
+ D. Greenman. 
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in the field on two occasions with success. The first was after the 1955 February 
1g earthquake near Quetta, when the mean epicentre and depth of focus of 
aftershocks were determined using records from the observatory and three field 
stations (Khan 1958). Of particular interest is the unusually shallow depth of focus 
(7km). This had been suspected to be true of earthquakes in this region on 
qualitative evidence (West 1935). The sedimentary rocks of the region are 
generally more than 7 km thick. 

The seismicity of Pakistan with special reference to the problem of designing 
seismic resistant structures has been reviewed (Khan 1956). 

At the request of Professor Antonio Marussi of Trieste, who visited the Insti- 
tute during the course of the Italian Expeditions in the Karakoram region, a 
review of the seismicity of the Hindukush and Pamir areas was made (Choudhury 
1959). ‘The Hindukush is an unusual deep focus (20c km) area, and the question 
was raised as to whether there is any form of coupling between this area and the 
adjacent normal focus area of the Pamirs. The epicentral map includes gravity 
anomalies known at the time. Some new work carried out with Professor 
Coloumb and Mme Labrouste in Paris is included in the review (Choudhury 
1956, 1958). Of special interest is his estimate of 75 km to the Mohorovitié dis- 
continuity in the Hindukush area using a reflected phase off the discontinuity 
from the deep focus shocks. 


Geomagnetism 


The observatory was installed at the Geophysical Institute in 1953 by Messrs 
S. Abu Tahir and S. A. A. Kazmi. The Chittagong equipment was installed for 
training, in a domestic house rented for the purpose, by Mr S. J. Ahmad, the 
present geomagnetician. A third station is planned for operation in the far north 
of West Pakistan, probably at Gilgit. 

At Quetta, base values were controlled by means of three QHMs and an Earth 
Inductor. Askania HTMs arrived last year to replace the QHMs. Ruska and 
Danish Quick-running magnetographs record the diurnal variations. Full details 
of the observatory, behaviour of the equipment, and the Annual Summary of 
observations up to the end of 1955 have been published (Kazmi 1956). The 
corpuscular eclipse, recorded at Quetta during the total sun eclipse of 1954 June 
30, has been described (Kazmi & Wienert 1955). An independent check on the 
precision of the Quetta Observatory up to 1955 has recently been made (Kazmi 
1959) in the course of a study of Sg variations under Dr Laursen at Copenhagen. 

The ionospheric equipment of the Geomagnetic Observatory is a standard 
NPL type which was installed by Messrs Ziauddin Ahmad and Musharraf Hussain. 
A three-mast rhombic aerial with vertical orientation is used. Geomagnetic and 
ionospheric bulletins are published monthly. 

The geomagnetic field survey of the country is being carried out by the Geo- 
detic Branch of the Survey of Pakistan in charge of Mr S. Q. Hasan. Forty new 
stations were observed in 1953 and 1954 and 3 stations of the former survey 
repeated by Mr Rajput. Some results are discussed later. 


Physical geodesy 


Mr. Hasan plans to have a gravity station at the corner of 10-mile squares 
throughout the country. Each station will be identified by a permanent plinth 
tied to the first-order level network. In 1955, 104 new stations in West Pakistan 
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and 67 in East Pakistan were measured with a Worden Geodetic gravity meter by 
Mr S. I. H. Musvi. The National Base value of gravity used in this work is that 
measured by Woollard and his colleagues who obtained independent values of 
978-9634 cm/s? and 978-9626 cm/s? at Karachi Civil Airport (Woollard et al. 1952). 
The mean value, 978-9630 cm/s”, has been adopted. It will be necessary to repeat 
all the old pendulum stations since some of those already checked are now beyond 
the desirable limits of precision (see for example Thirlaway 1959). A series of 
modern pendulum observations between Colombo and Gilgit is desirable for 
gravity meter calibration in the region. 


Both the field geomagnetic and gravity results are due to be published this 
year in the First Geodetic Report. 


Geodynamics 


Seismological, geomagnetic, gravity and geological observations in the Quetta 
re-entrant region are being integrated at the Geophysical Institute. 

Earthquake epicentres follow the strike of the Tertiary and the Secondary 
structures, the great majority being located south and east of a geanticline mapped 
by geologists of the Photographic Survey Corporation. The relatively few earth- 
quakes north and west of this structure are associated with a series of wrench faults 
trending NNE. (The largest of these, which can be traced for several hundred 
miles on air photographs, passes through the towns of Chaman and Nushki.) 
The formal projection of epicentres along any section normal to the re-entrant 
structure shows a concentration of seismicity south or east of the geanticline. 

The largest earthquake of recent years in the re-entrant is that of 1955 Feb- 
ruary 18, about 10 miles north of Quetta, at a depth of 7km (Khan 1958). The 
first motion recorded at Quetta was dilatation, as were those of all the Indian 
stations save Madras. Hamburg and Stuttgart recorded compression. Under- 
thrusting or sinistral wrench faulting near Quetta are consistent with the few 
observations that have been collected. 

A feature noted in the Quetta seismograms from some earthquakes within a 
few miles at Quetta is the unusually large vertical P amplitude. The simplest 
explanation is that the foci of these earthquakes are almost beneath Quetta. 
The S—P intervals will therefore measure approximately the depth of focus. This 
is usually less than 15km, in agreement with the findings of West and Khan 
already cited. 

Some 50 vertical force values of the Earth’s magnetic field were measured 
along the main Chaman-Sibi road (175 miles) through Quetta in 1953 and 1954 
by the Survey of Pakistan. There is a step of about 200y in the otherwise un- 
disturbed gradient of about 17y/mile throughout the traverse. This occurs 25 
miles south of Quetta at the head of the Bolan Pass leading down across the strike 
through the re-entrant structures. On the other hand, 12 vertical force values 
between Quetta and Kalat (80 miles) along the strike show no significant break 
in the regional gradient. 

Gravity anomalies observed by the Survey of India within the re-entrant struc- 
tures have the following features (see Survey of India Geodetic Reports for the 
gravity values). 

(a) The isostatic isoanomalies follow the geological strike of the mountain arc; 


(b) large negative isostatic anomalies (— 75 mgal) border the Indus valley boun- 
dary of the arc; 
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(c) the isostatic and free air anomalies increase positively with height (ranges: 
200ft to 800oft, —75 mgal to +45 mgal (isostatic), and —140mgal to 

+120 mgal (free air)). 

(d) the Bouguer anomalies are independent of height. (Mean value of about 

—150mgal over the same altitude range.) 

This behaviour is in marked contrast with that of the gravity field in other 
mountainous areas. In the easterly striking Mekran ranges of similar tectonic age 
to the south, for example, the expected behaviour of Bouguer and free air anoma- 
lies in a compensated region is observed (Thirlaway 1959). 

It is believed that the rapid increase in thickness of relatively light Quaternary 
and Tertiary rocks between the geanticlinal axis and the Indus Basin is the cause 
of these four unusual features of the gravity field. If the mean density of these rocks 
is of the order of 2-4 g/cm, and not 2-67 g/cm? as assumed in the isostatic reduc- 
tions, a variation in thickness up to a maximum of 30000 ft of low density rocks 
will remove most of the unusual features in the gravity field. An oil boring by 
Pakistan Petroleum Ltd., at Bannh (Sibi) near an isostatic anomaly of —75 mgal 
ended at a depth of 13001 ft in Siwalik shales. The figure of 30000 ft is greater 
than Glennie’s downwarps by a factor of two (Glennie 1956). However, general 
agreement has been reached on the existence of considerable crustal downwarps 
bordering the mountain arc of the Quetta re-entrant. Confirmatory seismic 
reflection work by oil companies is not yet available for release. 

The geology of the mountain arcs is everywhere exposed and the more diffi- 
cult to interpret for that reason. The oldest rocks exposed are Permo-Carboni- 
ferous along the line of the inferred geanticline. North and west of this line the 
sequence thickens rapidly as a shale facies; south and east as a limestone facies. 
Thrust plane traces are everywhere parallel to the sinuous strike of the arcs. 
Over-thrusting is everywhere towards the south and east. A complementary 
series of active wrench faults seem to be causing the deformation of the arc round 
Quetta. The larger set has been mentioned in connection with the seismicity of 
the area; recent movement is demonstrable. This sinistral series strikes NNE and 
the complementary dextral series WNW. A maximum principal stress directed 
horizontally in a direction S 20° E (or in the complementary direction) could pro- 
duce both sets of wrench faulting (Anderson 1951) and the effect of continued move- 
ment along both sets would be to deform a simple arc into a sinuous re-entrant. 

A model which fits the observed seismicity, gravity field, and geology has been 
constructed. This requires initial failure by thrusting (overthrust to south and 
east) in a manner identical to the classical island arc development, and subsequent 
deformation of the simple arc by wrench faulting. The required change of the 
intermediate principal stress to the vertical direction is brought about in a natural 
manner by the accumulation of mountain ranges during the initial thrust-type 
failure. The magnetic anomaly of 200 also fits this model if it is assumed (contrary 
to experience) that the thick sediments underlying the overthrust basement slice 
are slightly more magnetic than the underthrust basement. ‘The position, but not 
the sign, of the magnetic anomaly was predicted from the model. A means of 
testing the model is by fault displacement studies. 


Applied Geophysics 
Oil prospecting by private companies has resulted in large gravity coverage of 
the lower and middle Indus Basin and East Pakistan, including an off-shore survey 
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in the Bay of Bengal. Seismic reflection surveys are also extensive. Mr Menan 
Khan, Superintending Geophysicist of the Geological Survey, has undertaken 
the collection of gravity information released by the companies. 

Gravity and seismic surveys can assess groundwater resources in the western 
mountain valleys (Ahmad & Thirlaway 1958). It is often possible to distinguish 
between river and lake environment deposits by locating buried channels with 
gravity surveys. ‘The lake deposits are normally much less porous and therefore 
less productive of water. Seismic surveys are used to locate layering within the 
valleys. The structure of the valleys, bounded by limestone ranges, is rift-like 
with terminating cross faults. A few thousand feet of light material (presumably 
mainly Siwalik) has collected within the rifts; in the Quetta valley for example, 
there is some 4000 ft of relatively light material at its deepest area just north of 
the city. Maximum anomaly variations between the centre of the valleys and the 
denser rocky edges reach 20mgal with gradients of 1omgal/mile. ‘The density 
differential is about 0-3g/cm*. The gravity work has been extended to ccver 
several adjacent valleys (Mufti & Siddiqui 1959). 

Small haematite deposits have been located near Dalbandin with vertical mag- 
netic balances (Farah & Aslam 1958). Four thousand stations were observed on a 
200 ft grid. Shear zones were well delineated on this scale. The largest anomalies 
in the area were 7 000, and they appear from size, shape and drilling information 
to date, to be volcanic necks. ‘This was a combined operation by the Geological 
Survey and Geophysical Institute. 

Other work by the geophysical branch of the Geological Survey includes 
coalfield gravity and seismic surveys, a large number of electrical resistivity surveys 
for water, and dam site surveys. The Geophysical Institute and University of 
Lahore have been used on some of these operations. It is hoped to begin publica- 
tion of results in the near future. 

Dr. K. Helbig began systematic magnetic surveys with University of Panjab 
students over the Hidden Range this year. ‘This structure has already been well 
explored gravitationally and it is possible that oil company interests may again 
extend north to the Range. Other geophysical work in this area has been carried 
out by the PWD Groundwater Organization, a group formed a few years ago to 
study waterlogging and salinity problems caused by irrigation works. Very detailed 
electrical resistivity and well logging surveys were introduced by consultants of 


the International Co-operation Administration. The results are being prepared 
for publication. 


The Geophysical Institute, 
Quetta: 


1959 June. 
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Reports of Discussions 


Geophysical Discussions of the Royal Astronomical 
Society 


A discussion was held in the rooms of the Royal Astronomical Society on 1958 
November 28 at 16h 15m on Turbidity Currents: High-Speed Underwater 
Flow of Sediment Suspensions. It was under the Chairmanship of Dr G. E. R. 
Deacon, F.R.S. He opened the discussion by commenting on the prominence 
this problem of oceanography had acquired as a consequence of its being recognized 
that turbidity currents may have been responsible for the breaking of submarine 
communication cables. ‘The Chairman introduced each of the speakers in turn. 

Dr A. S. Laughton of the National Institute of Oceanography dwelt on three 
classes of evidence of turbidity currents from which our knowledge of them has 
derived. First, turbidity currents have been observed in lakes. In Lake Mead, 
formed behind the Hoover Dam on the Colorado River, currents have been obser- 
ved moving with speeds up to half a knot, one or two metres in thickness, and lasting 
several days. These currents carry a suspension of silt and clay which can reach a 
density of 2-05 and they originate in the turbid river which empties into the lake. 
Second, in the deep oceans turbidity currents are thought to be responsible for 
the deposition of relatively coarse (sand) layers of sediment far from land. Also, 
such currents appear responsible for the breaking of submarine cables. Heezen 
and Ewing calculated the speed of such a current from cable breaks following the 
Grand Banks earthquake of 1929. It appeared that a current with a speed of 60 
knots could have been responsible. From this study Kuenen calculated the size 
of the current and the quantity of material moved. Subsequent submarine coring 
in this area of the Atlantic by Ewing and Heezen confirmed the nature of the current 
through the nature of the sediment it transported. Similar cable breaks due to 
turbidity currents have been known to occur off the Algerian and Columbian 
coasts. Third, laboratory experiments can produce turbidity currents. It can be 
seen that the currents produce lateral as well as vertical variations in sedimentation. 
Kuenen, after many experiments in tanks and troughs, produced a formula, 
modified from the hydraulics for the flow of a river, relating the velocity to the 
depth, slope, and density of the current. The speaker ended by stressing the 
apparent importance of these currents in producing features of lake and ocean 
floors: the flatness of the bottoms and the movement there of coarse and poorly 
sorted materials. He called for more knowledge of their actions to better under- 
stand their mechanics. 

Dr R. Stoneley, F.R.S. of Cambridge University, spoke next. In discussing 
the suggestion that the breakages of submarine telegraph cables, following the 
Grand Banks earthquake of 1929 November 18, occurred through the passage of a 
turbidity current, he remarked that the hypothesis had the merit of being more 
free from ad hoc assumptions than some of the other explanations proposed. ‘The 
idea merited further examination from the hydrodynamical standpoint. In parti- 
cular, two aspects had to be considered, (a) whether the turbidity current would 
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spread longitudinally too greatly to be effective, (b) whether the motion of the 
current would persist over distances of the order of 1000km without being 
damped out by viscous forces. Under (a), the spread of the sediment-charged mass 
of water longitudinally had much in common with the problems of the surge 
following the breaking of a dam and of the subsidence of a volume of water, 
considered respectively by R. F. Dressler, and by W. G. Penney and C. K. 
Thornhill. A rough calculation based on the work of these authors suggested that 
the depth of the current as envisaged by Kuenen would be much too great; the 
hypothesis would be more workable if the current had a thickness of only a few 
metres. The question (b) depended on the viscous drag of a stream at high Rey- 
nolds numbers; the theoretical and experimental work of G. H. Keulegan on 
arrested saline wedges suggested a much lower drag coefficient than was appro- 
priate for lower speeds. Keulegan’s experiments showed that the logarithm of 
the drag coefficient depended linearly on the logarithm of the Reynolds number, 
and an extrapolation of this result indicated a thinner turbidity current, of vertical 
thickness about 5m. Such a thickness would help to meet the difficulty men- 
tioned under (a). Dr Stoneley said that further investigations of these hydrodyna- 
mical problems were very desirable, and examination was needed of the effective- 
ness of such a current in transporting and sorting sediments to the degree required. 

Dr T. H. Ellison of Manchester University was the final speaker in the formal 
portion of the discussion. He dealt with a mathematical consideration of the mixing 
of the turbulent flow with the clear water above thus diluting it and thickening the 
flow. Due to this thickening any filaments of the dense flow which entered the 
upper layer would be reincorporated and not escape from the current. The problem 
was to know when this dilution would halt the flow. An experiment by Turner 
consisting of injecting a two dimensional jet of fresh water into salt water had 
been analysed and it gave some estimates of the Richardson number. If this was 
low and combined with a high slope there would be much mixing. But with the 
small slopes occurring in the case of the Grand Banks turbidity current the Richard- 
son number was large. Using the figures of Kuenen the speaker showed the dis- 
tances and velocities of that current to be reasonable. The development presented 
had derived from a study of the migration of methane along the roofs of mine 
tunnels. 

In the general discussion which followed Dr K. E. Latimer made some obser- 
vations of concern to submarine cable engineers. A layer of sand one metre thick 
could be enough to bury a caple beyond the depth at which it might be recovered 
by a grapnel. A 3 mile length of cable with little slack might be broken by a 
current of 2-3 knots. Thus turbidity currents with velocities and densities as 
great as those mentioned in the discussion could not only break but sweep away 
the cable. Turbidity currents were also of interest in the matter of understanding 
the origin of submarine canyons. Comm. Cousteau had observed a turbidity 
current started accidentally. The material involved in the Grand Banks turbidity 
current had presumably come from the St. Lawrence River. This river provided 
an excellent example of entrained flow for sea water amounting to 100 times the 
river’s volume was carried 500 miles in a current of about 2 knots. Dr Latimer 
asked if such entrainment might explain sudden temperature variations observed 
by submarine cable engineers. He also asked if turbidity currents seemed to be 
acted on by the Coriolis force. 

Dr Strong described how in decanting, a more rapid separation of fluid and 
sediment could be effected by tipping the container to increase the surface and 











250 Geophysical discussion of the Royal Astronomical Society 


permit flow of the dense lower layer along the sloping side of the vessel. Such a 
continuous rather than an avalanche current was comparable to the difference 
between the currents discussed for the ocean floor and those formed where a 
river entered a sea. Snow avalanches had many similarities to the former type. 

In replying to Latimer Dr Laughton said that one difficulty in the problem of 
submarine canyons was to know what started them. A fault or fault zone seemed 
a possible explanation. Regarding the importance of the Coriolis force he observed 
that there was a tendency off the Pacific coast of North America for currents to 
turn left rather than right as they came out of submarine canyons. He said Mennard 
held the view that this was due to deposition causing a levee to be built up on the 
right-hand side. Latimer confirmed this observation. 

In a series of short comments concerning the breaking of submarine cables 
points stressed included that the tension in the cable was also partly due to the 
motion of the ship laying or grappling the cable and that a current acting over a 
long length of cable could produce a large tension. Laughton stressed that the 
central region of a turbidity current was one of erosion. Strong said that small 
movement in sediments could shear a steel tube. Mention was made of man-made 
currents following bomb explosions in a harbour. 

The Chairman asked Ellison if there was a layer of mixing in front of a snow 
avalanche. The reply was that while the speaker had not looked into the problem 
he wondered if the particles were not too dense for this to be so. 

The Chairman thanked the speakers and other contributors to the discussion 
and adjourned the meeting at 18h 05 m. 





A discussion was held in the rooms of the Royal Astronomical Society on 
1959 January 23 at 14h 15m on The Automatic Reduction of Geophysical 
Data, under the Chairmanship of Sir Edward Bullard, F.R.S. 

The Chairman began the discussion by answering the possible question of why 
one would wish for the automatic reduction of data, by saying that this was because 
of the size of the problems undertaken and the heavy burden of reducing some data 
by desk calculations. It gave greater accuracy and saved time for the scientist. 
He gave as an example his application of automatic reduction to magnetometer 
data from marine traverses obtained by Discovery IJ. ‘This method illustrated 
steps typical of many operations. In obtaining the data there was a chain consist- 
ing of the measuring equipment, an analogue-digital converter, a punch control 
and the punch. A chart recorder to plot the data for visual monitoring of the 
measurements was very helpful, as was the automatic input of additional digital 
data through the punch control. The speaker preferred to punch on to a tape 
because the order was kept clear. There was no shrinking to cause trouble in the 
automatic reading of the tape as was sometimes the case with cards and one num- 
ber was represented by one row of holes. In the automatic reduction process carried 
out on the EDSAC II machine, the date was first printed on the output, then the 
time was read, the magnetometer data read, the magnetic field calculated, head- 
ing and diurnal corrections applied, differences between the value determined and 














Geophysical discussion of the Royal Astronomical Society 251 


neighbouring values tested to reject erroneous observations, the ship’s position 
was interpolated and the output of the machine punched on tape. It was stressed 
that the punching of indicators on the tape could help by telling the machine what 
was to follow and by providing for contingencies. ‘The process of testing figures 
with some limits of variation was very useful for removing faulty observations and 
if these were printed separately, later inspection could ensure that valid data were 
not discarded. 

The Chairman then introduced the speakers and recommended that the discus- 
sion of each paper should immediately follow its delivery. 

Mr R. H. Merson of the Royal Aircraft Establishment spoke on the reduction 
of interferometer and kinetheodolite observations of Earth satellites. This work 
had been in progress since shortly after the launching of the first Russian satellite. 
Prediction of the position of the satellites had been of first concern but the study 
of their orbits was later extended to obtain information about the air density at 
high altitudes and about the gravitational field of the Earth. Various optical and 
radio methods of observation had been used; visual observation, radio Doppler 
studies, radio-interferometry, pilot balloon theodolite observations and kine- 
theodolite observations. ‘The last method gave an accuracy of +0-005 deg in 
angle measurements and of +0-o02s in time. The less accurate methods were still 
useful in obtaining initial estimates of the unknown elements of the satellite orbit. 
This orbit was then adjusted by successive approximations to give a least-squares 
fit with all the observational information. Automatic reduction of the data in 
a computer was necessary to achieve the speed needed in a prediction service. A 
Pegasus digital computer was used for the major part of the work. From an 
interferometer record a table of two shifted direction-cosines against time was 
obtained; the shifted direction-cosines were true direction-cosines plus unknown 
constants. These data were used to find the height, longitude and time as the 
satellite crossed the observer’s latitude, and the orbital inclination. In making 
this analysis the size and shape of the orbit (represented by the semi-major axis a 
and eccentricity e) were assumed to be known. a and e were obtained from the 
nodal period 7 and perigee distance r,,. 7 was provided by the prediction service 
and r,,, which changes comparatively slowly, from a combination of various types 
of observation. Initial estimates of the unknown elements were given to the 
computer which then worked out the values of the shifted direction-cosines in 
terms of these estimates, compared the results with the observed values, and used 
the residuals to improve the estimates by the method of differential correction. 
Kinetheodolite observations (elevation, azimuth, time) were converted to direc- 
tion-cosines and the analysis then proceeded along the same lines as indicated 
above. 

A later computer programme transformed the results into standard orbital 
elements. The behaviour of one of these, the Right Ascension of the ascending 
node §2,, had been used to estimate the values of coefficients in the expression for 
the Earth’s gravitational potential. From observations of Sputnik 2 (19578), 
combined with American results from their satellites 1958« and 195882, it had 
been deduced that 

J = (1625 + 0-001) x 1073; 
D = (8+3)x 10%. 


It followed that the Earth’s polar flattening e = 1/(298-15 + 0°09). 
During the discussion the Chairman asked if this slightly higher value of the 
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ellipticity over older results might be due to pear-shaped harmonics of the gravita- 
tional field or a poor distribution of the earlier results. Merson said he thought 
the latter seemed likely. Dr A. H. Cook suggested there might be poor statistical 
reduction of the older figures. 

The second speaker was Mr G. A. Corby of the Meteorological Office who 
discussed numerical weather forecasting with the aid of computers. The forecast- 
ing proceeded from the meteorological data recorded at various weather stations 
and compiled on maps. These data could have been selected from maps for 
computing by human agency but it had been thought more desirable to let the 
computing machine enter into even this part of the work. The machine took the 
data and interpolated between the observed points to obtain values of the variables 
at a gridwork of points evenly spaced over the area covered. A least squares reduc- 
tion was then begun to minimize the difference between observed and calculated 
values with a weighting of points to favour data from nearby meteorological stations. 
Some maps were shown which compared favourably the machine controlled maps 
with the hand-drawn ones. This project was still in the development stage. 

Questions were put to the speaker about transmission errors. He replied that 
there were two or three errors per 250 observations (the number of weather 
stations serving western Europe). Mr E. Gold said that when line telegraphy had 
been in use for transmitting observations, about one figure group in 2000 had 
been the frequency of transmission errors. This order of frequency for errors 
seemed inevitable and demanded a well-distributed network of observing stations 
for proper detection. He asked how long it took to prepare these charts with the 
aid of the computer. ‘The speaker replied that the time from beginning to end 
was one and a half hours. Dr Fox asked whether local irregularities in the observed 
data were not extremely troublesome as the numerical processes for solving 
equations appropriate to a selected model of the atmosphere could be somewhat 
unstable. Mr Corby replied that short-period disturbances could be troublesome 
to forecasting and a smoothed curve was most desirable. 

Dr M. H. P. Bott of Durham University spoke next about his use of a digital 
computer to reduce and interpret data about the gravitational field obtained by 
regional surveys over geological formations. One application had been to compute 
the correction due to the terrain of the area. Another was to calculate the gravita- 
tional field over irregularly-shaped bodies in the Earth’s crust. When these 
geological units could be divided into a set of prisms, or two-dimensional bodies 
whose boundaries were straight lines, the problem was attacked by relaxation 
methods. Assuming the density contrast between the body and the surrounding 
rocks to be known, the shape of the lower surface of the body could be calculated 
if the top surface was defined, by successive approximations. In the discussion 
which followed this paper the point was made that this knowledge about the 
density contrast was necessary if there were not to be multiple answers. 

Dr D. C. Martin, from the Royal Society, outlined the problems of data 
processing arising from the International Geophysical Year. Three regional 
World Data Centres were collecting the data from each subject of study. Some of 
this was to be published in the Annals of the I.G.Y., the rest would appear in 
established publications. Automatic treatment of the data was not planned for 
the whole of it but doubtless much would be handled this way by the investi- 
gators concerned with it. The Chairman pointed out that a general treatment of 
all the data by automatic methods had been discussed and a decision taken against 
it by the British National Committee for the 1.G.Y. 
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The fifth speaker that afternoon was Mr M. J. Tucker of the National Institute 
of Oceanography. His remarks referred to measuring and analysing the spectrum 
of random-gaussian processes of which sea waves and turbulence were typical 


examples. The basic limitation, independent of instrumentation, was expressed 
by the formula 


eBL = 1. 


Here, e was the proportional random error in the estimate of E(f), the power 
density of the spectrum; B was the bandwidth of the equivalent filter, the range of 
frequencies over which E(f) was smoothed and L was the length of the record or 
the time of observation. For a 10 per cent accuracy and 10 per cent frequency 
resolution, this required a record containing 1000 waves: a digital record might 
need to consist of 5000 numbers. In a naturally occurring phenomenon, the 
characteristics might have changed over a period of time of this length. By com- 
paring two related records (for example, surface wave height and the pressure on 
the sea bed immediately below), the random errors might be reduced or eliminated 
and higher accuracies be obtained. 

The speaker said analogue analysers were generally expensive to buy and 
cheap to run, whereas if a digital computer could be hired, digital methods involved 
little capital cost but a high cost per analysis. Digital methods were capable of 
greater accuracy. Digital methods that used correlation as the first step were 
economical of computer time and best for many purposes, though periodogram 
analysis and numerical filters could be useful in particular applications. 

In the discussion which followed this paper the point was raised that in some 
cases a natural phenomenon was stable enough for a sufficient length of time to 
justify the data being recorded as digits. The Chairman asked if it was not true 
that with digital Fourier analysis the time requirements would be similar to those 
of analogue methods. The speaker replied that the digital analysis required many 
more multiplication operations and was thus likely to need more time. Dr P. 
Vigoureux pointed out that for a disturbance which had appreciable periodic com- 
ponents a spectral representation could be obtained only by using a very narrow 
sampling band. In this case the measurement was no simpler than were auto- 
correlation methods. He stressed that the sampling band width should always be 
stated for such measurements. 

The possible use of a digital computer for ionospheric forecasting was the 
topic of a communication by Mr C. M. Minnis of the Radio Research Station. 
Such forecasting was concerned with making estimates of the highest radio fre- 
quencies which could be used in transmitting radio waves for long distances by 
reflecting the waves from the ionosphere. The first stage was to study records 
of certain ionospheric measurements and then, from a given trend of solar activity, 
predict future values of ionospheric parameters. The second stage was to compute 
the required radio frequencies which were functionally dependent on these para- 
meters. As an illustration of how a digital computer could help with these tasks, 
forecasting the critical frequency (f) of the F2 layer was described. Observation 
showed that the monthly median values of f were useful and that they depended 
on latitude and longitude as well as on the intensity of solar radiation, hour angle 
(H), and the mean solar declination (D). The solar intensity could be represented 
by a linear relation between the critical frequency and the ionospheric index 
(Ipe). By selecting values of H and D the solar intensity was expressed as a function 
of geographical location. Thus a large mass of hourly observations of f was 
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reduced to about 1 000 values (for a given month, from 65 locations) which could 
be dealt with by a computer at one time. The machine could then interpolate the 
functional dependence of f on Ip2 for any location and, given a forecast value of 
Ir2, compute a forecast value for f. In this way the derivation of the highest 
frequency which could be used for radio communication between two stations 
had been reduced to the evaluation of certain mathematical relations at a series 
of points along the great circle between these stations. 

Mr A. T. Dennison, British Petroleum Company, followed to speak about the 
use of digital computers in exploration geophysics. Applied to surveys of the 
gravity or magnetic fields, there were many operations of routine treatment of the 
data obtained by these surveys, similar to those described by Dr Bott, which could 
be treated by automatic reduction methods. Large storage capacity was perhaps 
the most restricting requirement for these computer applications. In seismic 
prospecting, a digital presentation of the data was not favoured. Here automatic 
treatment was usually concerned with corrections for one record at a time. It 
might be possible to compute automatically the correction required for the effect 
of the weathered layer on a seismic reflection from the records of first arrivals of 
the reflected waves. The corrected records were then combined for a visual 
presentation as profiles, variable intensity strips, or in some other manner. On 
the combined records a correlation from trace to trace was made to construct an 
interpretation of the subsurface structures. The judgment of which arrivals should 
be chosen for correlation was preferably left to an interpreter. Magnetic record- 
ing of the seismic traces was a relatively recent innovation to these methods. It 
permitted a greater range of frequencies to be recorded and in the subsequent 
automatic reduction of these data, various filters could be introduced. 

The eighth speaker was Dr G. W. Lennon of the Liverpool Observatory and 
Tidal Institute. He spoke on the processing of Earth tide data. This aimed at 
measuring the yield of the solid Earth to the gravitational attraction of the Sun 
and Moon and thereby obtaining information about the Earth’s elasticity. This 
study still required work on instrumentation as well as data-processing techniques. 
Records extending over long periods were difficult to obtain because of discon- 
tinuities. Many spurious contributions were difficult to separate from the tidal 
effects on the record. These could have a period very near that of the tide, as, 
for example, meteorological effects. An accuracy of 1 per cent was required in 
the basic values for the geophysical interpretation to proceed. In the work at 
the speaker’s Institute a modification of the methods developed by Dr Doodson 
was used to analyse 29 days’ observations. This was a manual technique involv- 
ing numerical filters to separate the tidal species, and correction of the individual 
constituents. Four reasons were given why digital computer techniques were 
unsuited to this work: (1) there was not a large amount of data but rather a collec- 
tion of small units; (2) the elimination of effects due to non-periodic movements 
had to be done at each stage of the analytical process; (3) in considering labour 
economy the analysis was the least demanding stage of the work; (4) it was 
desirable to preserve the non-periodic effects recorded for further geophysical 
investigations. 

The final speaker was Dr H. M. Iyer of India, temporarily at the National 
Institute of Oceanography. In the analysis of microseisms a three component 
seismograph station should provide information about the amplitude, period, 
spectra, direction of arrival, and composition of the impulses received. Dr Iyer 
described an automatic analyser which had been designed and applied to an 
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electronic feedback seismograph at the National Institute of Oceanography to do 
this analysis directly from the seismograph signals. The output of a seismograph 
was fed into a high torque galvanometer fitted with three systems of contacts. 
These contacts were used (1) to examine the position of the seismograph arm at 
short intervals, (2) to count the transits of the moving arm from the positive to the 
negative side during a known length of time, and (3) to examine whether two sets 
of waves were on the same, or different, sides of the zero line during a known 
length of time. This last technique was one due to Tomada and used for calculat- 
ing the correlation coefficients between the different seismic components in order 
to analyse the composition and direction of arrival of the microseismic waves. 
The speaker stressed that a treatment of the data such as this, which automated 
part of the analysis with simple equipment, could be of great help to scientists in 
countries such as his where the resources often did not permit more elaborate 
automation. 

The Chairman thanked the speakers and other contributors to the discussion 
and adjourned the meeting at 18h 20m. 





A discussion of papers received during the session was held in the rooms of the 
Royal Astronomical Society on 1959 March 20 at 16h 15m under the chair- 
manship of Dr R. Stoneley, F.R.S. 

The first speaker, Dr D. W. Allan, presented the paper by Mme H. A. 
Lubimova on The Thermal History of the Earth with Consideration of the Variable 
Thermal Conductivity of its Mantle (1958). Mme Lubimova had made an important 
contribution to studies of the Earth’s thermal history by making numerical calcula- 
tions which allowed for a thermal conductivity varying both with depth and with 
temperature (thus taking into account radiative transfer as well as normal thermal 
conduction). She found that the thermal history was still largely dependent on 
the initial temperature, the age of the Earth, and the distribution of radioactivity 
assumed, as with simpler mathematical models using constant conductivity. 
The course she proposed for the Earth’s thermal history was, very briefly: the 
Earth was formed by accretion of particles about 5 billion years ago. It was origi- 
nally cold but due to radioactive sources, kinetic energy of colliding particles and 
pressure, it heated up. Differentiation of radioactive materials was continuous due 
to partial melting but main differentiation occurred 1-6 billion years after its 
genesis. The increase of temperature was more rapid at greater depth than closer 
to the surface, but the core which was assumed as a mixture of iron and silicates 
started to melt only 1-2 x 10° years ago and is still melting. 

Due to a thermal conductivity minimum at a depth of 50-100 km heat generated 
in the centre could not escape. The depth of the region of escape was about 
sookm. The temperature inside this region had decreased about 100°C in the 
last billion years. The Earth as a whole could not melt in the process of its evolu- 
tion. 

Another possible thermal history for the Earth, which J. A. Jacobs and D. W. 
Allan discussed, assumed a slightly older Earth which was almost completely 
molten at 4°5x10° years ago causing differentiation of radioactive materials. 
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There was partial remelting near the surface later. The core was always molten 
and the inner core was formed at the beginning of the process of solidification. 

It was interesting that, on both these models, the Earth was heating up rather 
than cooling below 7oo or 800km depth. Might this be an explanation of the 
cessation of deep-focus earthquakes at that depth? 

Evaluating mathematical theories available at present, Dr Allan said that they 
had reached a necessary degree of complexity. However, more experimental work 
was needed, particularly on radioactive abundancies. More knowledge of the 
initial conditions of the solar system would be very helpful. It was imperative 
to decide how long ago the solar system was evolved and how long it took for 
the formation of the planets. In the discussion which followed Dr Runcorn said 
that he had been working recently on the mechanism of continental drift and 
that a continuous differentiation of radioactive material would help. For a model 
study, one must take into account convection currents in the core, and gravita- 
tion. Dr Benfield said that the surface of the Moon had not changed much and 
questioned whether it could be used to explain thermal history. Dr Allan 
answered that, according to his model, the Moon was solid throughout and cooled 
all the time. Dr A. H. Cook noticed that the curve of the melting point tempera- 
ture near the surface would allow for the presence of magma for volcanoes. Dr 
Stoneley mentioned that nobody attributed any anomalies on the seismograms 
to pockets of molten lava. 

The second paper presented was a joint paper by Mr L. R. Cooper and Mr C. 
Jones on The Determination of Virgin Strata Temperatures from Observations in 
Deep Survey Boreholes (1959). The knowledge of vertical temperature distribution 
in the ground around a favourable geological formation was of great importance in 
designing a new mine installation, a ventilation system and in evaluating future 
operational efficiency. The difficulty in obtaining temperature distribution was 
due to the fact that, once the ground had been disturbed, the virgin temperature 
changed. In order to measure the virgin temperature a special technique had to 
be developed including a specially designed thermal probe which could be used 
in conjunction with the drilling rig. The probe used a thermistor sensing device 
which was connected by a two-core cable to a Wheatstone bridge. The main 
advantage in using a thermistor probe was that a thermistor had an impedance of 
about 400kQ so that small variations in cable resistance due to changes in tem- 
perature along the borehole could be neglected. Some early measurements in 
coal measure rocks indicated that temperature geothermic gradient changes from 
100 to 40ft/°F. This was attributed to some low conductivity coal pieces in the 
strata pierced by the borehole. In the discussion, Dr T. F. Gaskell asked about the 
magnitude of the correction which had to be applied to measured temperatures in 
order to obtain true virgin temperature. Mr Cooper answered that the correction 
was 1 to 4 or 5°F, proportional to depth and depending upon the time necessary 
for withdrawal of the drilling rods. Dr W. H. Bullerwell mentioned that the Geo- 
logical Survey sometimes performed similar measurements but due to operational 
procedure it was not possible so far to pay much attention to the temperature of 
the bottom of the hole. Dr A. H. Cook suggested that the techniques described 
in this paper might be applied to heat flow measurements. 

The third paper of the session was presented by Dr A. H. Cook on Determina- 
tions of the Earth’s Gravitational Potential from Observations of Artificial Satel- 
lites (1958). If the Earth’s gravitational potential consisted of 1/r? terms only, the 
orbit of any artificial satellite would be a plane ellipse with focus in the centre of 
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the Earth. The fact that the satellite orbit rotated around the major axis was due 
to higher order terms of the potential. By observing the retrograde motion of the 
orbital plane the values of higher terms in the expression for potential could be 
evaluated. Utilizing observations of the Vanguard satellite, the following values 
had been obtained by Mr R. H. Merson: 


103] = 1°625+0-001 
10°9D = 8+3 
and for the ellipticity e~! = 298-2. 

Dr J. de Graaff-Hunter pointed out that the coefficients J and D, as found in 
geodetic work including gravity, were evaluated on the assumption that the level 
surface was a spheroid and D was chosen to fit the measurements. The spheroidal 
shape was worked out on the assumption that the Earth was a fluid rotating body 
and consequently we should not hang on to the old value of D when the new direct 
determinations of this factor gave a significantly different value. 

Dr R. L. G. Gilbert could not be present and his paper on An Investigation 
into the Calibration of Gravity Meters (1958) was presented by Dr W. H. Bullerwell. 
This paper described the work of Dr Gilbert at the Dominion Observatory, 
Ottawa. In large scale surveying with gravity meters it was very important to 
check if calibration constants of gravity meters changed. The repeated tests of 
two North American and one Worden gravity meter showed that the calibration 
constant changed over a period of time. It was concluded that measurements over 
the range of 1000 mgal might be preferable with pendulums. A calibration device 
built into the gravity meter which resembled the Askania gravity meter construc- 
tion did not operate as satisfactorily as was expected. Dr Bullerwell then described 
results of calibration tests done by the Geological Survey over a British calibration 
line. Six Worden instruments were tested. It was discovered that the change of 
calibration constant was proportional to the age of the instrument, increasing with 
time. This could be explained either by the ageing of the instrument or by a change 
in calibration procedure by the manufacturer. Dr Cook pointed out that an accur- 
acy in positioning the ball of 1 part in 108 was necessary in order to have a reliable 
built-in calibration device. 


The Chairman thanked the speakers and other contributors and adjourned the 
meeting at 18h 10m. 
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The Magneto-Ionic Theory and its Applications to the 
Ionosphere 


J. A. Ratcliffe 


(Cambridge University Press, 1959, 40s.) 


Mr J. A. Ratcliffe is famous, both nationally and internationally, for his capacity 
to select and simplify the fundamental features of complex phenomena, and his 
influence has stimulated the Cavendish Laboratory ionospheric group until it has 
become the Mecca of young research workers in this field. He has taken an active 
interest in the magneto-ionospheric theory since it was first developed and pub- 
lished one of the first detailed accounts of the results of the theory (Wireless 
Engineer 1933, 10, 354.) Even today, more than twenty-five years ijater, this 
paper can be recommended as an accurate and clear introduction to the subject. 
Throughout this period, Mr Ratcliffe has been particularly interested in the 
physical significance of the mathematics and has devoted much time and atten- 
tion to its demonstration. We may thus expect that his monograph on magneto- 
ionic theory will not only be a first-class contribution, but is likely to be remarkably 
clear and readable. Our expectations are not disappointed. 

The first impression given by this book is that it could well be taken as a model 
for similar works in other fields where progress is hampered by the difficulties and 
complexities of the basic mathematical theory. This is primarily due to the care 
with which the subject matter and degree of development have been limited so that 
the important features of the theory are brought out in simple form. At the same 
time, the intricacies of the theory are investigated sufficiently fully to remove the 
main difficulties of those who wish to proceed to a more exact and complete analysis 
of particular problems. The annotated bibliography is designed to guide this de- 
velopment and, in effect, adds several hundred pages of detailed analyses. Great 
care has been taken to underline the physical significance of the theory, the tech- 
nique of developing the theory from both macroscopic and microscopic points of 
view being particularly effective. However, possibly the most valuable feature 
of the book, both to undergraduate, young research worker, and expert is the 
collection and clarification of many difficulties and fallacies which have been en- 
countered in the past but are not readily available in the literature. 

The author’s experience of the difficulties met by research students in under- 
standing and using magneto-ionic theory is probably unequalled and has provoked 
him to develop new methods of explaining these difficulties. This is reflected in 
the book, which contains a large amount of valuable material which has never been 
published before. It is remarkable how little the standard figures and curves 
usually reproduced elsewhere contribute to the author’s argument. The biblio- 
graphy is arranged so that additional material can readily be found if needed. 

Roughly half the book is concerned with the development, properties and physi- 
cal significance of the magneto-ionic theory and roughly half with the problems 
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which arise when it is applied to the ionosphere. While, of course, much of this 
theory can be applied directly to non-terrestrial magneto-ionic media, such as the 
solar atmosphere, the development is primarily slanted to be directly applicable 
to the ionospheric problems of the interaction of a radio wave and the ionospheric 
media. In both cases the analysis corresponds to a ray treatment, the more com- 
plex phenomena associated with coupling, which need a full wave treatment, being 
discussed in outline but not in detail. References are given to the fuller treatments. 
With this limitation, which is imposed by the author’s decision to make his book 
simple and readable, the book is sound and complete, the only omission noted 
by the reviewer being the absence of any comment on, or reference to, Booker’s 
(J.1.E.E. 111 1947, 94, p- 171) analysis of propagation in terms of impedances. 
While, admittedly, this approach is fundamentally different to that selected by 
the author and has still to be expanded to be fully applicable to magneto-ionic 
topics, it provides yet another useful viewpoint the main points of which could 
have been summarized in an appendix. This would have been particularly valuable 
to radio propagation engineers, many of whom approach practical propagation 
problems in terms of telegraph theory. However, this attractive and detailed 
account of magneto-ionic phenomena in the ionosphere will undoubtedly remove 
many misconceptions. An excellent summary of the scope of the monograph can 
be obtained quickly from the list of contents though some readers may find certain 
of the last few sections to be slightly over-compressed. 

This is an exceptional book which can confidently be expected to become a 
standard text book on magneto-ionic theory as applied to the ionosphere. While, 
in the future, it may be desirable to add new sections, the main content of the 
book is not likely to date. There is no doubt that workers in this field will owe a 
considerable debt of gratitude to Mr Ratcliffe for this valuable monograph. 


W. R. Piccortrt. 


Fonctions Sphériques de Legendre et Fonctions Sphéroidales 
Louis Robin 


(Tome I 1957, 201 pp. + xxxv 4,000 FF. Tome II 1958, 384 pp. + viii 5,000 FF. 
Gauthier-Villars, Paris) 


The author of these two volumes, and a third which is not yet published, has 
set out to provide a comprehensive and up-to-date reference book on the proper- 
ties of Legendre Functions. The first volume begins with a useful detailed 
summary, under the heading “Introduction et vue d’ensemble”’, of the contents 
of all three volumes. 

Chapter I starts with the Legendre equation and the author proceeds to derive 
properties of the solutions P,(j) and On(u) for positive integral values of m and 
real or complex ». Chapter II extends this analysis to the associated Legendre 
equation and the associated Legendre functions P7'(«) and O”(y) for integral 
values of m and n, n being positive. 

Chapter III deals with integral spherical harmonics. Volume I ends with 


tables giving Pa(u), On(u), Pr(u) and O7(u) as functions of w for small integral 
values of m and n. 
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In Volume II the scope is greatly enlarged by extensions to complex values of 
mand n. In Chapter IV P™(u) and Q"(y) are defined by means of contour 
integrals for any real or complex m, n and uw. The properties of these generalized 
functions are then developed and exhibited as extensions of those for integral m, 
n. In Chapter V are obtained various asymptotic series for P™ and OQ” valid for 
large |m| or |n|. Chapter VI deals with the expression of a given function as a 
series in terms of associated Legendre functions. Volume III will deal with 
theorems of addition of Legendre functions, applications to ellipsoids of revolu- 
tion, Gegenbauer functions and spheroidal functions. 

These volumes are written with great skill and clarity; they form a compen- 
dium of information comparable with Watson’s Bessel Functions. In general full 
proofs are given, and the reader is assumed to be familiar with hypergeometric 
series, the method of steepest descents, asymptotic expansions and Lebesgue 
integrals. But for anyone needing only to extract a formula suitable for applica- 
tion in a problem of mathematical physics, the summary in the Introduction will 
form a convenient pointer, and the careful setting out of all formulae with stated 
conditions of validity should ensure correct use. The author is an engineer and 
the book arose from needs discovered in dealing with problems of fluid mechanics, 
geophysics and electromagnetism. Workers in such fields are likely to find it of 
very great value. ' 


E.R.L. 





Russian translations of summaries of original papers 


(Russian summaries of the foregoing papers are unavoidably held over until 
the next issue.) 


(MEPEBO)[bI PESIOME IIPEJIbIJYUIUX JOKIIAJOB HA PYCCKHA A3SbIK 
HEM3SBEsKHO 3AJEPKUBAIOTCA JO CIE YIOUJELO BBbIITYCKA,) 
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